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1.  INTRODUCTION 


In  this  report  we  assess  a  number  of  physical  and  chemical  processes 
for  their  effect  on  strlatlon  structuring,  with  particular  emphasis  on  phenomena 
present  at  times  after  a  nuclear  event,  but  prior  to  the  re-establishment  of 
the  undisturbed  ambient  magnetic  field.  During  this  time  Interval  the  plasma 
Is  far  from  thermal  equilibrium,  and  there  are  many  sources  of  free  energy  to 
drive  plasma  Instabilities.  One  could  attempt  a  detailed  plasma  physical 
description  of  the  phenomenology;  such  an  effort  even  If  theoretically  success¬ 
ful  would  be  difficult  to  verify  observatlonally. 

Alternately  we  have  sought  to  estimate  the  effect  at  time  t2  of  several 
physically  and  chemically  "well-established"  (e.g.,  non  plasma-wave)  mechanisms 
on  structure  of  specified  wavelength  and  amplitude  present  at  some  earlier 
time,  t^.1  As  time  progresses,  one  expects  the  plasma  free  energy  (and 
hence  the  variety  and  Intensity  of  Instability)  to  decrease.  It  Is  then 
possible  that  plasma  structure  at  t2  can  be  constructed  out  of  structure 
present  at  tj  which  evolves  by  non  plasma-wave  processes  between  tj  and  t2. 

If  plasma  wave  processes  are  completed  prior  to  tj  It  Is  likely  that,  at  most, 
they  affect  structural  dynamics  from  tj  to  t2  only  parametrically. 

This  approach  has  the  advantage  of  simplicity  and  of  certainty  In  the 
qualitative  aspects  of  Its  Individual  chemical  and  physical  components.  Upper 
bounds  on  each  Fourier  component  of  the  density  structure  at  t2  obtain  readily 
when  one  assumes  that  the  amplitude  of  the  component  at  t1  Is  of  order  the 
background  density  (l.e.,  relative  "perturbations"  of  order  unity). 

The  values  of  tj  and  t2  need  not  be  the  same  for  the  various  physical 
and  chemical  processes  Invoked  to  affect  structure.  In  Section  2  we  consider 
the  effect  of  charge  exchange  and  Impact  Ionization  on  structural  decay.  In 
this  case  tj  can  be  the  time  at  which  air  Ionization  pick-up  Is  completed  at 
a  location  at  the  burst  point  altitude  of  the  bomb,  and  t2  can  be  the  time  at 
which  the  resulting  energetic  air  Ions  have  completed  their  conjugate  region 
deposition.  (This  latter  time  Is  generally  prior  to  the  relaxation  of  the 
magnetic  field  to  Its  undl  turbed  st>  e.)  In  Section  3  we  consider  the  effect 
of  recombination  on  structu.«i  dec  /.  In  Section  4  we  consider  the  effect  of 


Coulomb  scattering  on  structural  decay.  For  these  two  cases  tj  can  be  the 
time  at  which  slow  (non  suprathermal  velocity)  ionization  is  created  and  t2 
can  be  the  time  of  re-establishment  of  the  undisturbed  magnetic  field. 

The  results  for  the  charge  exchange  and  impact  ionization  calculations 

I 

are  particularly  striking  in  that  slow  electron  structure  in  the  100-200  km 
altitude  range  due  to  fast  ions  created  at  an  altitude  of  400  km  exhibits  a 
sharp  decay  with  decreasing  wavelength  at  inverse  wavenumbers  below  20  kilo¬ 
meters  (see  Figures  11  and  13  for  details).  Such  information  provides  both 

2 

qualitative  and  quantitative  inputs,  "seeds,"  for  the  SCENARIO  code. 

In  addition  to  the  primary  results  of  Sections  2-4,  which  are  summar¬ 
ized  in  Section  5,  a  number  of  additional  results  relevant  to  striation  decay 
are  presented  in  Appendices  A-C.  Appendix  A  Is  a  theoretical  discussion  of 
the  relationship  between  magnetic  field  diffusion  and  particle  diffusion. 

This  provides  background  for  the  calculation  of  Section  4.  Appendix  B  is  a 
discussion  of  slow  ionization  production  by  12  keV  O”  ions,  as  occurring  in 
magnetic  storms.  Appendix  C  presents  new  results  on  the  suppression  of 
anomalous  diffusion  in  late  time  striation  structure  by  ionospheric  electric 
fields. 
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2.  THE  OECAY  OF  HIGH  ENERGY  IONIZED  STRUCTURE 
THROUGH  CHARGE  EXCHANGE  AND  IMPACT  IONIZATION 


The  undisturbed  Ionosphere  Is  a  complex  mixture  of  charged  and  neutral 
particle  species.  One  way  In  which  a  nuclear  detonation  can  modify  this  com¬ 
position  Is  through  the  generation  of  Ions  and  neutrals  with  energies  of  the 
order  of  10  or  even  100  keV.  Once  formed,  the  fast,  high  energy  neutrals  are 
transformed  Into  Ions  by  Impact  Ionization.  Between  charge  exchange  and 
Impact  Ionization  events,  high  energy  Ions  and  neutrals  can  Increase  the 
thermal  charged  particle  density  of  the  ionosphere  through  Ionization  colli¬ 
sions.  Eventually,  a  sufficiently  large  number  of  collisions  occur  so  that 
the  initially  high  energy  particles  become  part  of  the  thermal  background. 

The  purpose  of  this  section  is  to  quantify  the  degree  to  which  an 
initial  structure  consisting  of  high  energy  ions  is  degraded  as  it  flows 
down  field  lines  Into  denser  portions  of  the  ionosphere  and  to  determine 
the  density  of  thermal  Ionization  formed  as  a  result  of  ionization  by  high 
energy  ions  and  neutrals.  As  the  basis  for  determining  the  density  of  high 
energy  ions  and  neutrals  as  a  function  of  time,  the  following  equations  have 
been  solved: 


an, (r,t)  +  _ 

- - *  -VjUJn^r.t)  +  v2(t)n2(r,t)  (la) 

an2(r,t)  +  + 

- 3^ - *  -v2(t)n2(r,t)  +  v1(t)n1(r,t)  -  vQ(t)n2(r,t)  .  (lb) 


Here, 


n^(k,t)  »  high  energy  Ion  density 
n2(k,t)  *  high  energy  neutral  density 


vj(t)  »  rate  at  which  high  energy  Ions  are  destroyed  and  high 
energy  neutrals  are  created  as  a  function  of  time 


Vg(t)  =  rate  at  which  high  energy  neutrals  are  destroyed  and 
high  energy  Ions  are  created  as  a  function  of  time 

vQ(t)  *  rate  at  which  high  energy  neutral  structure  is  lost  by 
convection  or  diffusion. 


Writing  the  spatial  structure  of  nj(r,t)  and  n2(r,t)  in  terms  of 
Fourier  integrals. 


00 

n^r.t)  =  J  d3k  n1(r,t)  exp(ikx) 


[jt  +  vl^]  *  v2(t)n2(k,t) 

+  Mk,t)  =  vlni(^*t) 


(2a) 


n2(r,t)  =  j  d3k  n2(r,t)  exp(ikx)  (2b) 

•oo 

permits  Eqs.  (la)  and  (lb)  to  be  rewritten  in  the  following  way: 


(3a) 


(3b) 


In  writing  Eqs.  (1)  and  (3)  a  term  for  ions  comparable  to  vQ  for 
neutrals  has  not  been  included  since  structure  dimensions  larger  than  an  ion 
gyroradius  are  assumed.  Furthermore,  the  spatial  dependence  of  Vg,  Vj,  and 
v2  on  the  coordinate  orthogonal  to  the  field  lines  has  been  neglected.  This 
assumption  is  reasonable  if  the  spatial  variation  of  the  parameters  orthogonal 
to  the  field  line  is  neglected  and  the  evaluation  of  the  equations  is  in  a 
frame  of  reference  moving  along  the  field  lines  for  a  particular  velocity 
space  class  of  ions.  Hence,  the  solution  to  Eqs.  (1)  and  (3)  depends  upon 
the  Ion  velocity  considered. 

Two  forms  for  vQ  have  been  chosen.  For  the  case  when  the  neutral  mean 
free  path  Is  smaller  than  the  characteristic  dimensions  of  structure  a  diffu¬ 
sion  form  for  vg  Is  chosen.  In  Fourier  spacer^ 

v0  "  1(2  vni2/3v2  {4) 
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r 


where  vn^  is  the  component  of  neutral  velocity  perpendicular  to  the  magnetic 
field  lines.  For  neutral  mean  free  paths  larger  than  structure  dimensions  the 
diffusion  approximation  Is  not  valid  and  Vg  represents  a  convective  loss: 

Vg  -  k  vn/8  (5) 

where  v  Is  the  fast  neutral  speed  and  the  relation 
n 

Ax  Ak  =  \  (6) 

for  a  Gaussian  profile  has  been  used. 

In  addition  to  changing  form  as  either  fast  Ions  or  fast  neutrals, 
the  fast  particles  Ionize  the  cold  neutral  background  and  thereby  increase 
the  number  of  slow  ions  and  electrons.  The  expressions  giving  the  rate  of 
slow  electron  or  positive  Ionization  are 


dn  (k,t)  ,  _ 

— ^ —  =  n1(k,t)  ^  v$(t)n$(t)  +  n2(k,t)  £ 


vs'(t)ns(t) 


dn.  (k,t)  +  ___  +  _ 

— =  nj(k,t)  ^vs'(t)ns(t)  +  n2(k,t)  ^ 

s  s 


V(t)ns(t) 


(7a) 


(7b) 


where 

■  slow  electron  density 

■  slow  Ion  density 

■  collision  frequency  for  slow  electron  production  due  to 
fast  Ion  collisions  with  species  s 

v$'(t)  ■  collision  frequency  for  slow  Ion  production  due  to  fast 
neutral  collisions  with  species  s 


ne(k,t) 

n+(k,t) 

v$(t) 
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vs(t)  ■  collision  frequency  for  slow  positive  Ion  production  due 
to  fast  Ion  collisions  with  species  s 

v$'(t)  ■  collision  frequency  for  slow  positive  Ion  production  due 
to  fast  neutral  collisions  with  species  s 

In  writing  Eqs.  (7a)  and  (7b)  the  spatial  variation  orthogonal  to  the  magnetic 
field  of  vs(t),  v$'(t),  v$(t),  vs'(t).  and  n$  has  been  neglected. 

Eqs.  (3a)  and  (3b)  have  been  solved  with  k  assumed  to  be  one-dimensional 
and  just  after  a  nuclear  detonation  when  the  Ionosphere  Is  strongly  modified. 

MRC  code  results^  and  CIRA  tables^  have  been  used  as  the  basis  for  quantifying 
the  solutions  assuming  that  the  fast,  high  energy  species  Is  monatomic  oxygen. 
Also,  Reference  3  has  been  used  to  evaluate  the  collision  frequencies  and  vq. 

In  subsection  2.1  calculations  are  first  shown  for  the  400  km  burst 
altitude  Starflre  VER238  simulation  with  vertical  magnetic  field  lines,  and 
then  (assuming  negligible  change  In  the  Ionization  background)  for  the  same 
burst  with  magnetic  field  lines  making  an  angle  of  52.7°  with  the  vertical. 

In  subsection  2.2  calculations  are  shown  for  the  600  km  burst  altitude  TB600 
MRC  simulation.  In  subsection  2.3  calculations  are  shown  for  the  200  km 
altitude  LITTLE  ONE  simulation.  Comparison  of  the  results  of  the  various 
subsections  gives  Insight  Into  the  variables  affecting  Ionization  deposition 
and  structural  degradation. 

Magnetic  field  and  background  density  variation  used  In  subsection  2.1 
Is  shown  In  Figures  1  and  2.  The  variables  for  subsections  2.2  and  2.3  are 
shown  In  Figures  3  and  4,  and  5  and  6,  respectively.  Figures  1-6  are  micro¬ 
fiche  output  associated  with  Reference  5. 
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Figure  1.  Magnetic  field  geometry  for  the  400  km  altitude 
Starfish  VER238  simulation.  Fast  Ions  were 
assumed  to  originate  at  400  km  altitude  and  to 
move  down  along  a  magnetic  field  line  passing 
through  the  region  of  magnetic  field  compression. 
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Figure  3.  Magnetic  field  geometry  for  the  600  km  TB600 
MRC  simulation.  Fast  Ions  were  assumed  to 
originate  at  600  km  altitude  and  to  move  down 
along  a  magnetic  field  line  passing  through 
the  region  of  magnetic  field  compression. 
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Figure  6.  The  cold  neutral  background  appropriate  to  the 

200  km  Little  One  simulation.  For  this  simulation, 
with  (roughly)  constant  cold  neutral  mass  density 
at  altitude  below  200  km  along  a  magnetic  field 
line,  particles  were  assumed  to  move  along  the 
cold  neutral  mass  density  contour  for  3.0  x  10"15 
gm/cm3. 
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LITTLE  ONE  (200km) 


2.1  400  KM 

Based  on  the  MRC  simulation  of  STARFISH  at  400  km  and  the  geometry 
of  Figures  1  and  2,  calculations  have  been  carried  out  to  determine  the  Impact 
of  charge  exchange  and  Impact  Ionization  on  the  decay  of  fast  Ionized  struc¬ 
ture  moving  down  field  lines.  Also,  the  slow  electron  and  Ion  generation 
resulting  from  Ionization  by  fast  Ions  and  neutrals  has  been  evaluated.  The 
only  differences  with  the  calculations  of  Section  2.2,  q.v.,  are  that  the 
calculations  for  STARFISH  were  started  at  400  km  and  the  CIRA  atmosphere 
chosen  was  appropriate  to  STARFISH.  The  results  are  tabulated  In  Tables  1 
and  2  and  Figures  7  to  14.  Table  1  shows  the  slow  Ionization  density  pro¬ 
duced  by  fast  particles  derived  from  the  original  fast  Ion  for  an  Inverse 
wave  number  of  500  km.  Results  are  presented  for  electrons  and  Ion  produced 
by  160  keV  and  30  keV  Ions  and  neutrals.  Table  2  shows  the  field  line  Inte¬ 
grated  production  of  slow  electrons  and  Ions  as  a  function  of  Inverse  wave- 
number  for  both  160  keV  and  30  keV  cases. 

Figures  7  to  9  show  fast  Ion  and  fast  neutral  density  as  functions 
of  wavenumber  and  altitude  (or  time).  Figure  10  shows  the  Inverse  wavenumber 
variation  of  fast  Ions  and  neutrals  at  an  altitude  of  110  km.  Figures  11  and 
12  show  respectively  the  slow  electron  and  slow  positive  Ion  production  at 
specified  wavenumbers  relative  to  the  production  for  Inverse  wavenumber  of 
500  km  by  fast  particles  at  160  keV.  Hence  the  curves  are  Indications  of 
the  relative  values  of  Fourier  components  of  density  structure  If  one  assumes 
a  500  km  scale  length  background.  For  a  background  of  smaller  scale  length 
the  ordinate  for  the  same  abscissa  would  be  greater.  Figures  13  and  14  are 
similar  to  11  and  12  but  for  fast  particles  of  30  keV. 

The  qualitative  results  are  similar  to  those  of  Section  2.2;  however, 
because  of  the  different  cold  neutral  atmospheres,  there  Is  a  greater  degree 
of  slow  Ionization  at  relatively  short  wavelengths  (e.g.,  5  km)  produced  by 
the  STARFISH  than  for  the  600  km  burst  (Tables  2  and  4).  Note  also  that 
slow  electron  and  Ion  production  for  short  wavelengths  Is  greater  for  160 
keV  particles  than  for  30  keV  particles. 

The  MIC  simulation  assumes  that  the  magnetic  field  Is  aligned  with 
gravity.  However,  the  magnetic  field  lines  In  the  real  Ionosphere  actually 
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Table  1.  The  slow  electron  production  by  160  keV  particles,  n  *®®(k_), 

160  e  ® 

the  slow  Ion  production  by  160  keV  particles,  n.  (kA),  the 

30  ” 

slow  electron  production  by  30  keV  particles,  ne  (^0),  and 

the  slow  Ion  production  by  30  keV  particles,  n.30(k„),  are 

'  2i 

tabulated  as  a  function  of  altitude.  kQ  *  1/500  km  .  All 
quantities  are  In  units  of  particles/cm  per  single  fast  Ion 
at  400  km  altitude. 


H(km) 

"e160(ko)xl°6 

n+ieo(k0! x 106 

ne30(ko)xl06 

n+30(ko)xl06 

400 

3.5 

X 

10"  3 

3.6 

X 

10“2 

1.5 

X 

10"  3 

2.9 

X 

10"3 

390 

7.0 

X 

10“  3 

4.6 

X 

10"2 

3.2 

X 

10"3 

4.3 

X 

10"3 

380 

1.1 

X 

10"2 

5.3 

X 

10"2 

5.3 

X 

10"3 

5.4 

X 

10“3 

370 

1.9 

X 

CM 

1 

o 

H 

6.6 

X 

10"  2 

9.0 

X 

10"3 

7.7 

X 

10‘3 

360 

2.9 

X 

10"Z 

7.7 

X 

10"2 

1.3 

X 

10"2 

9.2 

X 

10"3 

350 

4.6 

X 

10“2 

9.5 

X 

10"2 

2.1 

X 

10"2 

1.2 

X 

10"2 

340 

6.5 

X 

10"2 

1.1 

X 

10"1 

3.0 

X 

10"2 

1.4 

X 

10"2 

330 

9.7 

X 

10"2 

1.4 

X 

10"1 

4.4 

X 

10"2 

1.8 

X 

10"2 

320 

1.3 

X 

10'1 

1.6 

X 

10"1 

5.8 

X 

10‘2 

2.0 

X 

10‘2 

310 

1.9 

X 

10"1 

2.0 

X 

10"1 

8.1 

X 

10"2 

2.6 

X 

10"2 

300 

2.5 

X 

10"1 

2.5 

X 

10" 1 

1.0 

X 

10"1 

3.0 

X 

10"2 

290 

3.4 

X 

10" 1 

3.2 

X 

10“ 1 

1.4 

X 

10" 1 

4.1 

X 

10"2 

280 

4.6 

X 

10" 1 

4.3 

X 

10"1 

1.8 

X 

10"1 

5.8 

X 

10"2 

270 

6.3 

X 

10"1 

5.9 

X 

10" 1 

2.3 

X 

10"1 

8.5 

X 

10"2 

260 

8.8 

X 

10"1 

8.0 

X 

10"1 

3.1 

X 

10"1 

1.3 

X 

10"1 

250 

1.2 

X 

10° 

1.1 

X 

10° 

4.1 

X 

10"1 

1.9 

X 

10" 1 

240 

1.7 

X 

10° 

1.6 

X 

10° 

5.6 

X 

10"1 

2.9 

X 

10"1 

230 

2.5 

X 

10° 

2.2 

X 

10° 

7.6 

X 

10" 1 

4.3 

X 

10"1 

220 

3.6 

X 

10° 

3.2 

X 

10° 

1.0 

X 

10° 

6.6 

X 

10"1 

210 

5.3 

X 

10° 

4.6 

X 

10° 

1.4 

X 

10° 

1.0 

X 

10° 

200 

7.8 

X 

10° 

6.8 

X 

10° 

2.1 

X 

10° 

1.6 

X 

10° 

190 

1.2 

X 

101 

1.0 

X 

101 

3.0 

X 

10° 

2.4 

X 

10° 

180 

1.8 

X 

101 

1.6 

X 

101 

4.4 

X 

10° 

3.9 

X 

10° 

170 

3.0 

X 

101 

2.5 

X 

101 

6.8 

X 

10° 

6.4 

X 

10° 

160 

4.9 

X 

101 

4.2 

X 

101 

1.1 

X 

101 

1.1 

X 

101 

150 

8.6 

X 

101 

7.3 

X 

101 

1.8 

X 

101 

1.9 

X 

101 

140 

1.6 

X 

102 

1.4 

X 

102 

3.4 

X 

101 

3.7 

X 

10l 

130 

3.6 

X 

102 

3.1 

X 

102 

7.3 

X 

101 

8.2 

X 

101 

120 

1.0 

X 

103 

8.4 

X 

102 

1.8 

X 

102 

2.3 

X 

102 

20 


Figure  7.  For  160  keV  particles,  the  fast  Ion  density  as  a  function  of 
wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density  as 
a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized  to 
the  Initial  fast  Ion  density  at  time,  t*0  sec,  and  altitude  of 
400  km,  and  plotted  versus  time,  t(sec),  and  altitude,  H(km). 
The  Initial  neutral  density  Is  zero.  The  different  types  of 
lines  used  In  the  figure  represent  different  wavelengths: 

k”*  ■  5  km - ,  ■  20  km - ,  ■  100  km  —  • — ,  and 

■  500  km . . 
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Figure  8.  For  30  keV  particles,  the  fast  ion  density  as  a  function  of 
wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density  as 
a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized  to 
the  initial  fast  ion  density  at  time,  t  =  0  sec,  and  altitude  of 
400  km,  and  plotted  versus  time,  t(sec),  and  altitude,  H(km). 

The  initial  neutral  density  is  zero.  The  different  types  of 
lines  used  in  the  figure  represent  different  wavelengths: 


n,  k,t=0 


110  H(kro) 


Figure  9.  For  30  keV  particles,  the  fast  Ion  density  as  a  function  of 
wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density  as 
a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized  to 
the  initial  fast  Ion  density  at  time,  t  ■  0  sec,  and  altitude  of 
400  km,  and  plotted  versus  time,  t(sec),  and  altitude,  H(km). 

The  Initial  neutral  density  Is  zero.  The  different  types  of 
lines  used  In  the  figure  represent  different  wavelengths: 


1  fast  Ion  and  fast  neutral  density  at 
km  normalized  to  the  Initial  fast  ion 
versus  Inverse  wavenumber  in  kilometers 


Figure  11.  For  160  keV  particles,  the  slow  electron  production  as  a 

function  of  wavenumber  and  time  Is  normalized  to  the  corres¬ 
ponding  quantity  In  Table  1  (kQ"*  *  500  km)  and  plotted  versus 
altitude,  H(km),  and  time,  t(sec).  The  different  types  of 
lines  in  the  figure  represent  different  wavenumbers:  k'*  ■ 


Figure  12.  For  160  keV  particles,  the  slow  positive  Ion  production  as  a 
function  of  wavenumber  and  time  Is  normalized  to  the  corres¬ 
ponding  quantity  in  Table  1  (k^1  ■  500  km)  and  plotted  versus 
altitude,  H(km),  and  time,  t(sec).  The  different  types  of 
lines  In  the  figure  represent  different  wavenumbers:  k"*  ■ 


Figure  13.  For  30  keV  particles,  the  fast  electron  production  as  a 

function  of  wavenumber  and  time  Is  normalized  to  the  corres¬ 
ponding  quantity  In  Table  1  (k  «  500  km)  and  plotted  versus 
altitude,  H(km),  and  time,  t(sec).  The  different  types  of 
lines  In  the  figure  represent  different  wavenumbers:  k’*  » 

5  km - .  ■  20  km  -  -  ■ 


•Ti, 


»♦  t“o-t) 


.6  t(sec) 
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Figure  14.  For  30  keV  particles,  the  fast  Ion  production  as  a  function 
of  wavenumber  and  time  Is  normalized  to  the  corresponding 
quantity  In  Table  1  (kQ“*  *  500  km)  and  plotted  versus 
altitude,  H(km),  and  time,  t(sec).  The  different  types 
of  lines  In  the  figure  represent  different  wavenumbers: 
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make  an  oblique  angle  with  gravity.  To  determine  the  effect  of  magnetic  field 
geometry  on  slow  electron  and  Ion  production,  the  relative  production  rates  of 
slow  electrons  and  Ions  for  magnetic  field  lines  which  make  a  realistic  angle 
of  52.7°  with  gravity  have  been  calculated.  The  production  rates  for  slow  Ions 
and  electrons  for  this  angle  have  been  plotted  In  Figures  15  to  17  for  three 
altitudes  (l.e.,  250,  200,  and  150  km,  respectively)  relative  to  the  production 
rates  for  magnetic  field  lines  which  are  aligned  with  gravity.  The  figures 
Indicate  that  the  slow  particle  production  rates  are  decreased  for  the  oblique 
angles,  especially  for  larger  wavelengths  and  lower  altitudes  primarily  because 
with  oblique  angles  there  Is  more  time  for  the  neutral  convective  or  diffusive 
loss  processes  to  decrease  the  number  of  particles  available  for  slow  particle 
production. 
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Figure  15. 


The  rates  of  slow  electron  and  Ion  production  by  160  and  30 
keV  particles  for  magnetic  field  lines  at  an  angle  of  52.7° 
with  the  vertical  relative  to  production  rates  for  magnetic 
field  lines  at  an  angle  of  0°  with  the  vertical  [i.e..  Re*®® 
ne160  (e  =  52.7)/ne160  (e  *  0),  Re30  3  ne30  (e  «  52.7)/ne30 
(e  3  0)]  are  plotted  versus  Inverse  wavenumber  at  an  altitude 
of  250  km.  The  Initial  Ion  and  neutral  densities  at  400  km 
altitude  are  taken  to  be  1  and  0,  respectively. 
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Figure  16.  The  rates  of  slow  electron  and  Ion  production  by  160  and  30 
keV  particles  for  magnetic  field  lines  at  an  angle  of  52.7° 
with  the  vertical  relative  to  production  rates  for  magnetic 


field  lines  at  an  angle  of  0°  with  the  vertical  [i.e.,  ReADU  = 
ne160  (e  «  52.7)/ne160  (9*0),  Re30  =  ne30  (e  =  BZ^J/ne30 
(9  *  0)]  are  plotted  versus  Inverse  wavenumber  at  an  altitude 
of  200  km.  The  Initial  Ion  and  neutral  densities  at  400  km 
altitude  are  taken  to  be  1  and  0,  respectively. 


Figure  17.  The  rates  of  slow  electron  and  ion  production  by  160  and  30 
keV  particles  for  magnetic  field  lines  at  an  angle  of  52.7° 
with  the  vertical  relative  to  production  rates  for  magnetic 
field  lines  at  an  angle  of  0°  with  the  vertical  [i.e..  Re160 
ne160  (8  =  52.7)/ne160  (e  *  0),  Re30  =  ne30  (e  *  52.7)/ne30 
(e  =  0)]  are  plotted  versus  Inverse  wavenumber  at  an  altitude 
of  150  km.  The  Initial  Ion  and  neutral  densities  at  400  km 
altitude  are  taken  to  be  1  and  0,  respectively. 


2.2  600  KM 

In  Figures  18-20  and  Tables  3  and  4,  the  results  of  our  calculations 
for  fast  particle  structure  degradation  and  slow  particle  production  are 
graphed  and  tabulated  for  Initial  fast  Ion  energies  of  160  keV  arid  30  keV. 

The  Initial  fast  Ion  density  at  an  altitude  of  600  km  at  time  equal  to  zero 
Is  normalized  to  be  one  for  all  wavenumbers  while  the  Initial  fast  neutral 
density  Is  set  equal  to  zero.  The  fast  Ion  velocity  makes  an  angle  of  60° 
with  the  magnetic  field.  As  the  fast  Ions  move  down  the  field  line  they  are 
converted  Into  fast  neutrals  by  charge  exchange.  The  neutrals  subsequently 
either  again  revert  to  fast  Ion  form  due  to  Impact  Ionization  or  are  lost  as 
a  result  of  convection  or  diffusion. 

The  following  figures  and  tables  specifically  address  the  600  km 

5 

VER248  MRC  simulation.  Based  on  the  simulation,  the  magnetic  field  Is 
assumed  to  be  vertical,  aligned  with  the  gravity  vector,  while  the  cold 
neutral  background  Is  assumed  to  be  undisturbed  by  the  nuclear  burst  to 
lowest  order.  These  variables  are  shown  In  more  detail  in  Figures  3  and  4. 

In  Figures  18  and  19  normalized  n..(k,t)  and  nn(k,t)  are  plotted  versus  alti¬ 
tude  and  time  for  160  keV  particles.  The  curves  demonstrate  that  the  fast 
ion  contribution  to  structure  decays  most  strongly  for  short  wavelengths. 

The  neutral  density  starts  off  at  an  Initial  value  of  zero  and  reaches  a 
maximum  for  short  wavelengths  and  then  decays.  For  longer  wavelengths  the 
neutral  density  does  not  reach  a  maximum  but  continually  Increases  for  the 
times  and  altitudes  of  Interest.  Figures  20  and  22  are  comparable  to 
Figures  18  and  19  except  that  the  particle  energy  Is  30  keV.  The  general 
qualitative  features  ascribed  to  the  earlier  figures  also  apply  to  Figures 
20  to  22.  In  Figures  23  and  24  the  total  fast  Ion  and  neutral  density  are 
plotted,  at  310  km  and  120  km,  respectively.  It  Is  evident  that  a  greater 
number  of  fast  Ion  and  neutral  particles  are  lost  from  the  short  wavelength 
as  compared  to  the  long  wave  contribution  to  structure.  This  result  Is  im¬ 
plied  by  the  neutral  diffusion  and  convection  models  which  are  most  effective 
for  short  wavelengths.  It  Is  also  noteworthy  that  for  long  wavelengths  when 
the  convective  and  diffusive  losses  are  weakest  the  density  profiles  evolve 
In  a  quasi-steady  state  manner  with 

vjn^k.t)  ■  v2n2(k,t)  .  (8) 
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Table  3.  The  slow  electron  density  production  by  160  keV  particles, 

k_ ) ,  the  slow  Ion  density  production  by  160  keV 
“  o  160 

particles,  n  (k  ),  the  slow  electron  production  by 
0  30 

30  keV  particles,  ng  (k  ),  and  the  slow  Ion  production 
by  30  keV  particles,  n+^(kQ),  are  tabulated  as  a  function 
of  altitude.  kQ  =  1/500  km.  Quantities  are  in  units  of 
particles  produced/cm  per  single  particle  at  600  km  altitude. 
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X 

10-2 

3.7  x  10'2 

420-410 

3.0  x  10’1 

2.8  x  10’1 

1.1 

X 

10'1 

4.2  x  10“2 

410-400 

3.5  x  10'1 

3.2  x  10’1 

1.2 

X 

10"1 

5.0  x  i»)”2 

400-390 

4.0  x  10_1 

3.7  x  10"1 

1.4 

X 

10'1 

5.8  x  10'2 

390-380 

4.6  x  10'1 

4.2  x  10’1 

1.6 

X 

10" 1 

7.0  x  10*2 

380-370 

5.3  x  10’1 

4.8  x  10’1 

1.7 

X 

10’1 

8.1  x  10“2 

370-360 

6.1  x  10_1 

5.5  x  10'1 

2.0 

X 

10_1 

9.7  x  10~2 

360-350 

7.0  x  10-1 

6.3  x  10’1 

2.2 

X 

10"1 

1.1  x  10'1 
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Table  3  (Continued) 


H(kra)  ne160(ko)  xlO6  n+160(kQ)  xlO6  ne30(ko)  x  106  n+30(ko)xl06 

350-340  8.2  x  10"1  7.3  x  10"1  2.5  x  10'1  1.4  x  10’1 

340-330  9.4  x  10"1  8.4  x  10*1  2.9  x  10'1  1.6  x  10"1 

330-320  1.1  x  10°  9.9  x  10"1  3.3  x  10'1  1.9  x  10"1 

320-310  1.3  x  10°  1.1  x  10°  3.7  x  10"1  2.3  x  10'1 

310-300  1.8  x  10°  1.6  x  10°  4.9  x  10'1  3.2  x  10'1 

300-290  2.1  x  10°  1.8  x  10°  5.7  x  10'1  3.9  x  10’1 

290-280  2.5  x  10°  2.2  x  10°  6.6  x  10"1  4.7  x  10_1 

280-270  3.0  x  10°  2.6  x  10°  7.8  x  10"1  5.7  x  10_1 

270-260  3.6  x  10°  3.1  x  10°  9.2  x  10'1  7.0  x  10"1 

260-250  4.3  x  10°  3.8  x  10°  1.1  x  10°  8.6  x  10_1 

250-240  5.3  x  10°  4.6  x  10°  1.3  x  10°  1.1  x  10° 

240-230  6.5  x  10°  5.7  x  10°  1.6  x  10°  1.3  x  10° 

230-220  8.2  x  10°  7.1  x  10°  1.9  x  10°  1.7  x  10° 

220-210  1.0  x  101  9.0  x  10°  2.4  x  10°  2.2  x  10° 

210-200  1.4  x  101  1.2  x  101  3.1  x  10°  2.9  x  10° 

200-190  1.8  x  101  1.6  x  101  4.0  x  10°  3.9  x  10° 

190-180  2.5  x  101  2.1  x  101  5.4  x  10°  5.4  x  10° 

180-170  3.5  x  101  3.0  x  101  7.6  x  10°  7.7  x  10° 

170-160  5.3  x  101  4.5  x  101  1.1  x  101  1.2  x  101 

160-150  8.6  x  101  7.3  x  101  1.8  x  101  1.9  x  101 

150-140  1.5  x  102  1.3  x  102  3.1  x  101  3.4  x  101 

140-130  3.2  x  102  2.7  x  102  6.4  x  101  7.3  x  101 

130-120  1.0  x  103  8.4  x  102  1.9  x  102  2.3  x  102 
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Table  4.  The  total  production  of  slow  electrons  and  ions  between 
600  and  120  km  per  single  fast  ion  at  600  km  altitude 
is  tabulated  as  a  function  of  inverse  wavenumber. 
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Figure  18.  For  160  keV  particles,  the  fast  Ion  density,  as  a  function 

of  wavenumber  and  time,  r^(k,t),  and  the  fast  neutral  density 
as  a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized 
to  the  Initial  fast  ion  density  at  time,  t*0  sec,  and  altitude 
600  km,  and  plotted  versus  time,  t(sec),  and  titude,  H(km). 
The  Initial  neutral  density  is  zero.  The  different  types  of 
lines  used  in  the  figure  represent  different  wavelengths: 


Figure  19.  For  160  keV  particles,  the  fast  ion  density,  as  a  function 

of  wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density 
as  a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized 
to  the  initial  fast  Ion  density  at  time,  t  =  1.368  sec,  and 
altitude,  310  km,  and  plotted  versus  time,  t(sec),  and  altitude 
H(km).  The  Initial  neutral  density  is  zero.  The  different 
types  of  lines  used  In  the  figure  represent  different  wave¬ 
lengths:  k"1  »  5  km -  ,»  20  km  -  -  -  ,  «  100  km  —  •  —  • 

and  «  500  km  . 
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Figure  20.  For  30  keV  particles,  the  fast  ion  density,  as  a  function 

of  wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density 
as  a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized 
to  the  initial  fast  ion  density  at  time,  t  *  0  sec,  and  altitude 
600  km,  and  plotted  versus  time,  t(sec),  and  altitude,  H(km). 

The  Initial  neutral  density  Is  zero.  The  different  types  of 
lines  used  in  the  figure  represent  different  wavelengths: 
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Figure  21.  For  30  keV  particles,  the  fast  Ion  density,  as  a  function 

of  wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density 
as  a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized 
to  the  Initial  fast  Ion  density  at  time,  t  =  0  sec,  and  altitude 
600  km,  and  plotted  versus  time,  t(sec),  and  altitude,  H(km). 

The  initial  neutral  density  Is  zero.  The  different  types  of 
lines  used  In  the  figure  represent  different  wavelengths: 
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Figure  22.  For  30  keV  particles,  the  fast  Ion  density,  as  a  function 

of  wavenumber  and  time  n^(k,t),  and  the  fast  neutral  density 
as  a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized 
to  the  initial  fast  Ion  density  at  time,  t  *  1.368  sec,  and 
altitude,  310  km,  and  plotted  versus  time,  t(sec),  and  altitude 
H(km).  The  Initial  neutral  density  is  zero.  The  different 
types  of  lines  used  In  the  figure  represent  different  wave¬ 
lengths:  k“*  ■  5  km - ,  *  20  km  -  — ,  ■  100  km  —  •  —  •, 

and  ■  500  km  . 


*e  23.  The  total  fast  Ion  and  fast  neutral  density  at 

H  =  3i0  tan  normalized  to  the  Initial  fast  Ion  density 
versus  Inverse  wavenumber  In  kilometers. 
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The  slow  electron  and  Ion  production  due  to  Ionization  by  160  keV 
and  30  keV  Ions  and  neutral  particles  Is  tabulated  In  Table  3  for  the  500 
km  wavelength  and  plotted  In  Figures  25  -  28  for  shorter  wavelengths.  In 
the  figures  the  slow  electron  and  Ion  densities  are  normalized  to  the  corres¬ 
ponding  500  km  wavelength  values  found  In  Table  3.  Because  the  calculations 
are  based  on  the  CIRA  tables  which  are  tabulated  for  10  km  distance  Intervals, 
the  numbers  In  Table  3  are  actually  the  number  of  particles  which  would  be 
generated  for  a  time  Interval  corresponding  to  the  traversal  of  160  or  30  keV 
particles  across  ten  kilometers  of  altitude.  To  determine  an  average  value 
for  slow  electron  and  Ion  density  per  cubic  centimeter  over  a  ten  kilometer 
distance,  the  numbers  in  the  table  are  divided  by  10®.  There  are  four  signi¬ 
ficant  qualitative  features  In  Table  3  and  Figures  25  -  28.  First,  the  number 
of  slow  particles  produced  Increases  at  lower  altitudes  for  the  500  km  wave¬ 
length  calculation  primarily  because  the  neutral  convection  or  diffusion 
losses  are  relatively  weak  compared  to  shorter  wavelengths  and  because  the 
cold  neutral  background  is  greater  at  lower  rather  than  at  higher  altitudes. 
Second,  relative  to  the  500  kilometer  slow  electron  and  ion  generation,  the 
generation  at  shorter  wavelengths  becomes  increasingly  weaker  especially  at 
lower  altitudes.  The  reason  for  this  effect  is  the  loss  of  fast  ions  and 
neutrals  due  to  the  convective  and  diffusive  neutral  loss  processes  at  short 
wavelengths.  Third,  the  total  density  of  slow  electrons  and  ions  which  would 
be  contained  within  one  cubic  centimeter  If  the  slow  particle  generation  had 
taken  place  within  one  cubic  centimeter  can  be  calculated  by  adding  up  the 
contributions  In  each  column  of  Table  3.  It  is  clear  that  carrying  out  the 
calculation  down  to  an  altitude  of  310  km  still  permits  the  creation  of  a 
larger  number  of  electrons  and  ions  than  were  In  the  original  500  km  struc¬ 
ture  at  an  altitude  of  600  km  (see  Table  4).  However,  for  shorter  wavelengths 
the  number  of  slow  electrons  and  Ions  is  smaller  than  for  k“*  =  500  km.  This 
effect  is  again  due  to  the  loss  mechanism  which  especially  affects  short  wave¬ 
length  structure.  Fourth,  the  160  keV  particles  generate  a  greater  number  of 
slow  electrons  and  ions  than  do  the  30  keV  particles. 
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Figure  25.  For  160  keV  particles,  the  slow  electron  production  as  a 

function  of  wavenumber  and  time  Is  normalized  to  the  corres¬ 
ponding  quantity  In  Table  3  (k0_1  *  500  km)  and  plotted  versus 
altitude,  H(km),  and  time,  t(sec).  The  different  types  of 
lines  In  the  figure  represent  different  wavenumbers:  k-1  a 


Figure  26.  For  160  keV  particles,  the  slow  positive  ion  production  as  a 
function  of  wavenumber  and  time  Is  normalized  to  the  corres¬ 
ponding  quantity  In  Table  3  (kQ-1  =  500  km)  and  plotted  versus 
altitude,  H{km),  and  time,  t(sec).  The  different  types  of 
lines  In  the  figure  represent  different  wavenumbers:  k"1  * 


Figure  27.  For  30  keV  particles,  the  fast  electron  production  as  a 

function  of  wavenumber  and  time  Is  normalized  to  the  corres¬ 
ponding  quantity  In  Table  3  (kQ-*  *  500  km)  and  plotted  versus 
altitude,  H(km),  and  time,  t(sec).  The  different  types  of 
lines  In  the  figure  represent  different  wavenumbers:  k"1  ■ 


Figure  28.  For  30  keV  particles,  the  fast  positive  ion  production  as  a 
function  of  wavenumber  and  time  Is  normalized  to  the  corres¬ 
ponding  quantity  In  Table  3  (ko“*  *  500  km)  and  plotted  versus 
altitude,  H(km),  and  time,  t(sec).  The  different  types  of 
lines  in  the  figure  represent  different  wavenumbers:  k"1  ■ 

5  km - ,  ■  20  km  -  -  *  100  km  —  • —  • 


2.3 


200  KM 


Calculations  to  determine  the  decay  of  fast  Ionized  structure  down 
field  lines  as  a  result  of  charge  exchange  and  impact  ionization  have  been 
carried  out  for  a  burst  altitude  of  200  km  with  the  ambient  atmosphere  chosen 

5 

from  Reference  6,  consistent  with  MRC  simulations.  Also,  based  upon  the  MRC 
simulations,  LOW  BEAR  200  and  LITTLE  ONE  200,  the  magnetic  field  lines  In  the 
vicinity  of  the  burst  are  taken  to  be  at  45  degrees  with  the  vertical  and  the 
cold  neutral  density  along  the  magnetic  field  lines  is  assumed  to  be  a  con- 
stant  and  appropriate  to  the  bubble  surface  (i.e. ,  3  x  10  grams/cm  ). 
Results  are  summarized  in  Table  5  and  Figures  29  to  35.  The  geometry  for 
this  case  Is  shown  in  Figures  5  and  6.  Because  contours  of  constant  density 
are  roughly  parallel  to  the  magnetic  field  and  the  cold  neutral  density  varia¬ 
tion  across  the  magnetic  field  was  not  taken  into  account  in  the  calculations, 
the  results  presented  here  actually  represent  an  underestimate  of  the  actual 
Impact  of  charge  exchange  and  impact  ionization  on  the  decay  of  the  initial 
fast  ion  structure. 

Neutral  effects  appear  singularly  weak  in  this  case  (compared  to 
those  previous)  because  the  high  degree  of  ionization  associated  with  a 
fireball  has  eroded  the  cold  neutral  density  appropriate  to  the  ambient 
undisturbed  ionosphere.  In  reality,  far  from  the  burst  point  the  cold 
neutral  density  approaches  ambient  and  slow  electron  and  ion  production 
becomes  much  larger  than  the  initial  fast  ion  density,  and  is  comparable 
to  the  lower  altitude  production  for  the  600  km  and  400  km  bursts  presented 
above.  This  is  not  determined  quantitatively  from  the  microfiche  data  of 
Reference  5. 
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Figure  29.  For  160  keV  particles,  the  fast  Ion  density  as  a  function  of 
wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density  as 
a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized  to 
the  initial  fast  Ion  density  at  time,  t  *  0  sec,  and  altitude, 
200  km,  and  plotted  versus  time,  t(sec),  and  altitude,  H(km). 
The  initial  neutral  density  Is  zero.  The  different  types 
of  lines  In  the  figure  represent  different  wavelengths: 

k“*  *  5  km -  ,  *  20  km  -  -  -  ,  3  100  km  —  •  —  • ,  and 

»  500  km  .  .  The  total  cold  neutral  density  is  taken 

to  be  3  x  10" 15  g/cm3. 
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Figure  30.  For  30  keV  particles,  the  fast  ion  density  as  a  function  of 
wavenumber  and  time,  n^(k,t),  and  the  fast  neutral  density  as 
a  function  of  wavenumber  and  time,  nn(k,t),  are  normalized  to 
the  initial  fast  ion  density  at  time,  t  *  0  sec,  and  altitude 
200  km,  and  plotted  versus  time,  t(sec),  and  altitude,  H(km). 
The  initial  neutral  density  is  zero.  The  different  types 


of  lines  in  the  figure  represent  different  wavelengths: 


K 


M 


The  total  fast  Ion  and  fast  neutral  density  at 
H  ■  100  km  normalized  to  the  Initial  fast  Ion 
density  versus  wavenumber  In  kilometers. 


Figure  32.  For  160  keV  particles,  the  slow  electron  production  as  a 
function  of  wavenumber  and  time  normalized  to  the  Initial 
fast  Ion  density.  The  different  types  of  lines  represent 

different  wavenumbers:  k“*  *  5  km  -  ,  =  20  km  -  -  -  , 

■  100  km  —  •  —  • .  and  *  500  km .  . 


Figure  33.  For  160  keV  particles,  the  slow  Ion  production  as  a 

function  of  wavenumber  and  time  normalized  to  the  Initial 
fast  Ion  density.  The  different  types  of  lines  represent 

different  wavenumbers:  k"*  ■  5  km -  ,  *  20  km  -  -  - 

■  100  km  —  •  —  •«  and  ■  500  km .  . 


Figure  34.  For  30  keV  particles,  the  fast  electron  production  as  a 
function  of  wavenumber  and  time  normalized  to  the  Initial 
fast  Ion  density.  The  different  types  of  lines  represent 

different  wavenumbers:  k"*  »  5  km -  ,  ■  20  km  -  -  - 

■  100  km  —  •  —  •,  and  ■  500  km .  . 
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Figure  35.  For  30  keV  particles,  the  fast  Ion  production  as  a 

function  of  wavenumber  and  time  normalized  to  the  Initial 
fast  Ion  density.  The  different  types  of  lines  represent 

different  wavenumbers:  k_1  »  5  km -  ,  »  20  km  -  -  - 

■  100  km  —  •  —  • , 


*  500  km 


3.  EFFECT  OF  OECAY  BY  RECOMBINATION  ON  THE  PSD 


All  Ionization  decay  appears  to  be  nonlinear  In  the  Ionization  density, 
since  at  some  point  In  the  recombination  process  (which  can  be  succession  of 
Interactions)  two  charged  particles  produce  a  neutral.  This  can  be  seen  In 
detail  In  References  7  and  8. 

Ionization  decay  should  obey  an  equation  of  the  form 


dn 

Ht  " 


(9) 


with 


a  -  decay  coefficient 
n  «  electron  density. 

We  will  do  specific  calculations  with  a  *  constant,  although  generally 
^  >  0.®  Usage  of  a  ■  constant  In  place  of  ^  >  0  would  tend  to  present 
a  lower  bound  on  the  smoothing  effect  of  nonlinear  decay,  but  further 
Intricacy  Is  possible. 

Consider  variation  over  a  time  Interval  with 


0  <  t  <  t  , 

(10a) 

n(t  *  0)  =  nQ  , 

(10b) 

n(t)  »  n 

(10c) 

Then  Integrating  Eq.  (9)  results  in 


or 


n  1  +  at  n, 


(lib) 


On  Fourier  transforming  both  sides  over  a  periodicity  length  X: 

-X 


n(k)  •/  e-”*’  .i.  dr 


1  +  at  nr 


(12) 


For  nn  at  <  1, 


n(k)  =  nQ(k) 


(13) 


Clearly,  nothing  much  happens  chemically. 
For 


nQ  at  »  1 


(14) 


one  has: 


Then: 


1  *  °*  "o  '  "o  0  ♦  S^t> 


«*  [*  '  V*  +  (n0lt)2  +  "•] 


(15) 


n(k  ■  0)  -  ^ 


X 

at 


(16a) 


but: 


n(k  f  0) 


(at) 


A 

l  f  _i_ 

J  nn 


e"1kX'  dX' 


(16b) 
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Combination  of  Eqs.  (16a)  and  (16b)  results  In  the  Interesting  conclusion 
that  the  ratio  n(k  f  0)/n(k  ■  0)  decreases  with  Increasing  time.  To  see 
this,  consider  the  case  (kj  f  0): 


n0  *  nQ(k  a  0)  +  2nQ(k1)  costkjX) 


,  k,  *  0 


2nQ(k1)  «  nQ(k  -  0) 


(17) 


Then: 


1  1  2n0(k.) 

=  1 - -  .  0  1  ..  cos(k  x) 

n0  nn(k  =  0)  n  Z(k  -  0)  1 


(18) 


and,  using  (16b): 

n(kx)  = 


1 


/ 


,  ,  2"0<kl>  1kl*' 

o  /  (-)  — j - e 

(atr  J0  n0Z(k  =  0) 


cos  kj X *  dX' 


X  n0^kl^ 


(at)2  n^(k  =  0) 


(19) 


Next  using  Eqs.  (16a)  and  (lib): 


n(kj)  ! 


nn(ki ) 


nn(*ci ) 


O'  1'  "0','1/  "0 

"^rr^  ‘  „  2<k  -  0)  ‘  "otk  •  ®>  "o<k  ■  °>  • 


(20) 


Associating  ng(k  =  0)  with  ng  results  In 


n^kl^ _ _no(kl>  noiklJ 


n(k  -  0)  "  nQ  nQ(k  «  0)  <<:  nQ(k  -  0) 


(21) 


Alternatively,  again  associating  nQ(k  *  0)  with  nQ;  n(k  ■  0)  with  n: 

_2  n0(k. ) 

ntfc.)  .  J-  -S-5-  . 

1  nQ  nQ 


(22) 
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Combination  of  Eqs.  (lib),  (22),  and  (13)  leads  to  Figure  36.  It 
Is  Interesting  to  note  that  n(kj)  Is  bounded  by  nQ*(*  1/at)  for  all  nQ; 
further,  that  If  nQ  >  nQ*.  n(kj)  varies  as  nQ*  (n0*/ng),  l.e. ,  n(kj)  decays 
as  t"2,  even  though  n  decays  as  t”1. 

Let  us  assume  that  Initially  (at  times  less  than  1  second  after  burst) 
one  has  structure  with  n0(kj)/nQ  -  1,  and  that  the  electron  temperature  varia¬ 
tion  among  bursts  Is  less  significant  than  the  density  variation.  Then  at 
times  after  burst  of  the  order  of  5  seconds,  one  would  expect  density  fluc¬ 
tuations  to  be  less  pronounced  at  densities  with 


n0 


> 


1 

5T 


-12  3  11  -3 

With  a  =  10  cm  /sec,  this  condition  becomes  ng  >  2  x  10  cm 


(23) 


Figure  36. 
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Assuming  full  Ionization  and  an  electron  temperature  of  five  times  the 

12  -3 

ambient,  this  corresponds  to  an  ambient  background  density  of  10  cm  , 
l.e. ,  significant  structural  non-sharpness  at  burst  altitudes  below  about 
100  km.  Data  from  Bluglll  and  Klngflsh  exhibit  structural  non-sharpness, 
while  Checkmate  and  Starfish  do  not  appear  to  exhibit  structural  non-sharp¬ 
ness.  This  is  in  keeping  with  our  theoretical  expectation. 

The  applicability  of  nonlinear  chemical  damping  of  structure  as  a 

source  of  structural  diffusiveness  is  confined  to  altitudes  below  110  km  if 

one  is  limited  to  times  of  the  order  of  5  seconds  or  less.  However,  this 

limitation  may  well  be  arbitrary.  One  can  also  consider  the  fate  of  structure 

generated  by  some  source  of  instability  (which  is  "turned  off"  at  time  tj)  in 

the  time  interval  between  tj  and  tg,  with  t2  denoting  the  time  for  onset  of 

structural  development  appropriate  to  the  gradient  drift  instability,  for 

instance.  If  one  takes  t0  -  300  seconds  and  t.  *  1  second,  the  density  limi- 

^  9  -3  ^ 

tations  apply  to  nQ  as  low  as  3  x  10  cm  .  With  temperatures  of  0.5  eV,  this 
would  affect  structures  up  to  an  altitude  of  160  km;  at  lower  temperatures, 
which  could  be  more  appropriate  at  times  of  the  order  of  hundreds  of  seconds 

Q 

after  burst,  from  Sappenfield's  a-coeffi dents  one  sees  that  the  effect 
could  be  even  more  extensive  in  altitude. 

One  can  calculate  the  effect  of  the  recombination  on  a  general  profile 
nQ(X)  through  Eq.  (12).  Roughly  speaking,  for  nQ(X)  <  1/at  there  is  no  loss 
of  structure;  for  nQ  >  1/at  the  structure  is  lost  more  rapidly  than  the  mean 
density.  To  get  some  feeling  for  the  effect,  we  consider 

n0  *  n0(k  =  [l  +  (1  "  £)  cos  kix]  »  0  <  e  <  1  .  (24) 


Then  for  nQ  at  »  1  from  Eq.  (16b),  one  has 


n(mk^)  = 


Ut)7 


n0U  -  6) 


j  -imkjX' 

(1  +  (1  -  e)  cos  kjX']  e 

m  *  integer 


(25) 
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On  substituting  z  *  k^X',  and  using  some  trigonometry  and  a  standard  formula: 

.2m  „ 


10 


•  Epiffinsr 


(-)Vm 


1 


(at)  /2  -e  i /c 


(26) 


with 


r  = 


/2  -  e  -  /e 


•/2TT~rD 


(27) 


One  can  verify 


0  <  r  <  1 


for 


0  <  e  <  1  . 


From  Eq.  (26),  an  exponential  spectrum  results  from  the  profile  of 
Eq.  (24).  The  variation  as  e  -*■  0  is  not  strongly  divergent.  Hence,  even 
for  rather  large  initial  modulations,  structure  at  kj  Is  lost  more  rapidly 
than  the  mean  density,  while  at  the  same  time  bounded  higher  harmonic  behavior 
Is  created. 

From  this  analysis  one  concludes: 


1. 


2. 


Chemical  recombination  produces  a  spectrum  with 


"0tfc  *  °> 

"ou  *  01  n0 


for  ng  at  »  1.  Put  another  way,  for  ng  at  »  1,  n(k  f  0)  -  t 
while  n(k  *  0)  ~  t”^. 


From  Figure  36,  for  all  possible  nQ  with  nQ(k  f  0)/nQ(k  ■  0)  fixed 
and  of  order  unity,  the  maximum  n(k  f  0)  at  time  t  is  of  order  nQ*, 
(*  1/at)  and  derives  from  nQ  *  nQ*.  For  nQ  >  n0*,  n(k  /  0)  < 
no*(no*/n0). 
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3. 


For  rig  at  »  1  and 

n0  =  n0(k  =  °)  u  +  (1  -  e)  cos  kjX]  ,  0  <  e  <  1  , 

an  exponentially  decaying  (in  wavenumber)  spectrum  of  higher  harmonics 
is  generated  [see  Eq.  (26)]. 
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4.  COULOMB  DIFFUSION  OF  FAST  IONS  AND  ELECTRONS 


In  this  section  expressions  for  the  Coulomb  diffusion  of  suprathermal 
electrons  and  ions  across  a  slab  geometry  magnetic  field  are  derived.  The 
plasma  model  used  for  the  calculation  Is  Identical  to  the  one  used  in  Appendix 
A  with  the  exception  that  an  x-directed  ambipolar  electric  field  will  be  re¬ 
quired  for  the  present  discussion.  Only  elastic  collisions  are  considered. 

The  equations  appropriate  to  the  description  of  transport  of  fast 

*•  i  H-13 

particles  are: 


j.  x 
e  o 


en  E  +  v  j‘ 
e  e  e  e 


(28a) 


Ji  x  B  -  m. 

0  .  - - VP.  +  en.E  -  Vi  —  ^ 


(28b) 


In  Eqs.  (28a)  and  (28b)  the  notation  is  similar  to  Appendix  A  except 
that  quantities  here  refer  to  fast  ions  and  electrons.  Also,  in  the  equations 
and  v.j  are  the  collision  frequencies  for  momentum  loss  with  the  less  ener¬ 
getic  background  particles  of  the  energetic  electrons  and  ions,  respectively; 
however,  the  reaction  between  energetic  particles  is  neglected. 

To  lowest  order  in  the  assumed  small  parameters  ve/|fle|  and 
the  solutions  to  Eqs.  (28a)  and  (28b)  are: 


je(0>  ■  y 


I—  — 

U  3X 


+  en  E  •  x 
e 


(29a) 


j  (°)  =  y  /  C  _L 

Jl  y  ax 


en^E  •  x 


(29b) 


To  next  order  in  the  small  parameter, 


V"  ■  ^  fe] 

1  (tT  *  enef  •  «) 

(30a) 

t  (1)  .  -  vimi  / c 

Ji  «  w 

)  (w  -  enfE  •  '*)  • 

(30  b) 

Because  of  quasineutrality. 

ne~  n. 

(31) 

it  is  necessary  that  there  be  no 

x-directed  current: 

0  =  ( je  +  jj)  *  x 

• 

(32) 

Equation  (32)  implies  that 

*  vimi(3Pi/3x) 

-  VgfaPg/ax) 

(33) 

eng(v 

eme  + 

and 

!e(1)  • ;  -  -  V11 

.  ViVt  j 

veme  +  vimi  e 

(34) 

Several  features  of  Eqs. 

(33)  and  (34)  are  noteworthy. 

First,  it  is 

possible  to  determine  a  constraint  on  the  ambipolar  electric  field  in  Eq.  (33) 
which  is  consistent  with  quasineutrality  if  we  scale  the  electric  field  as 


*  T 

E  *  x  ~± 


(35) 
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where  T  is  a  temperature  and  L  Is  a  scale  length.  From  Gauss'  law  this  con¬ 
straint  Is: 


T  ,  (V)‘ 

4„„A2  VL/ 


«  1 


(36) 


In  Eq.  (36),  the  parameter  aq  Is  the  Debye  length.  Second,  Eq.  (34)  is  inde> 
pendent  of  the  Ion  collision  frequency  if 


veme  <<:  vimi  *  (37) 

Physically,  this  corresponds  to  the  ion  transport  across  the  magnetic  field 
being  constrained  by  the  slower  electron  transport  rate.  Contrarily,  if 

vimi  <<e  veme  *  (38) 

the  electron  transport  is  constrained  by  the  slower  ion  transport  rate. 

To  complete  the  description  of  electron  and  ion  transport  across  the 
magnetic  field,  a  precise  form  for  the  electron  and  ion  collision  frequencies 
for  suprathe’Tnal  electrons  and  ions  is  necessary.  Such  expressions  are  found 
in  Ref.  14  under  the  assumption  that  the  thermal  particles  are  electrically 
charged.  For  fast  electrons  collisions  with  both  background  electrons  and 
ions  contribute  and 

vg  3  <7-7  \,e  +  3-9  xei}  x  10"6  ne  e-3/2  .  (39) 

In  Eq.  (39)  the  parameters  n  and  e  are  the  background  electron  density  (in 
-3  ® 

units  of  cm  )  and  energy  (in  units  of  electron  volts),  respectively.  The 
Coulomb  logarithms  Agg  and  Ag^  are  a  function  of  the  background  electron 
temperature  (Tg),  background  ion  temperature  (T^),  background  electron  den¬ 
sity,  and  the  relative  ion  to  hydrogen  mass  ratio,  p: 

A  ■  23  -  In  (ne1/2T /3/2)  T  <  10  eV 

ee  e  e  ,  e 

-  24  -  In  (n  1/2T  "X)  ,  T.  >  10  eV  ,  (40a) 

6  6  ® 
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Afii  *  23  -  In  (ne1/2Te"3/2)  ,  10  eV  >  Tg  >  T^/m, 

*  24  -  In  (ne1/2Te_1)  ,  Tg  >  10  eV  >  T^/n^ 

*  30  -  In  (n11/2T1“3/V1)  ,  Ti  >  T^/n^  .  (40b) 

The  ion  collision  frequency  depends  on  whether  the  energetic  ions  have  speeds 
which  are  in  excess  of  or  less  than  the  electron  thermal  speed.  For  ions 
with  speeds  which  are  larger  than  the  electron  thermal  speed  electrons  are 
the  primary  cause  of  drag  on  energetic  ions  and 

=  1.7  x  10‘4  y'1/2  Afii  n1  e'3/2  .  (41) 


where  u  is  the  ratio  of  ion  the  proton  mass.  However,  if  the  electron  thermal 
speed  is  larger  than  the  energetic  ion  speed  drag  on  both  electrons  and  ions 
contribute  to  the  slowing  down  of  a  fast  ion,  i.e.. 


^  =  1.8  x  10“7  W‘1/Z  Xi1ni  e"3/2  +  1.6  x  10"9  y_1  n£ 


Aei  Te 


■3/2 


(42) 


In  Eq.  (42) 


Ai i  =  41.5  -  0.5  In  (ne/Te)  (43) 

and  It  has  been  assumed  the  ion-ion  collisions  are  betweer  like  ion  species. 

On  using  Eqs.  (40a),  (40b)  and  (43)  to  compare  (39)  and  (41)  or  (42), 
one  readily  establishes  that  to  within  factors  of  unity 

>  meve 

Hence  particle  fluxes  can  be  estimated  through  the  simplified  form  of  Eq. 
(34): 
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Vex  ’  Vlx  -  -7^"  lx  (Pe  *  P1>  ’ 


(44) 


If  temperature  variation  Is  neglected,  the  density  diffusion  coefficient, 
Dn  Is  given  by: 


me 


Dn  "  “ITT  ve  (Te  +  V  * 
e  Bo 


(45) 


Let  us  seek  the  bound  on  wavelengths  such  that  there  Is  a  density 
decay  of  at  least  e~*  In  5  minutes.  This  condition  Is: 

D„  k  2  >  3  x  10"3  sec-1 
n  y  — 


or  equivalently  (for  Tg  =  T^): 


x  < 


2-5  x  10*&f  te)  • 


(46) 


It  Is  clear  that  significant  Coulomb  diffusion  perpendicular  to  B,  as  an 
early  time  phenomenon  (<.  5  minutes  time  scale)  Is  limited  to  wavelengths 
of  order  several  kilometers  or  less. 
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5.  SUMMARY  AND  OVERVIEW 


In  time  sequence  one  expects  first  charge  exchange  and  Impact  Ioniza¬ 
tion  (Section  2),  and  then  recombination  (Section  3),  and  Coulomb  diffusion 
(Section  4)  to  affect  the  Fourier  amplitudes  In  wavenumber  of  the  charged 
particle  density.  The  relationship  between  Inverse  wavenumber  (k“*)  and 
structure  scale  length  perpendicular  to  B  (Ln)  is  ambiguous  because  Ln  Is 
not  precisely  defined;  typically  however  one  expects  for  reasonable  defini¬ 
tions  of  L  that 
n 

L  -  V?  If1 
n 

(To  see  this  consider  a  wave  acos  kx,  a  «  1,  on  a  background  of  size 
unity.  Let  us  find  the  scale  length  corresponding  to  this  wave. 

To  first  order  In  o  we  have 

vn/n  «  k  acos  kx 

or  for  the  root  mean  squared  value 

(Vn/nl  m  c  =  (/dx  k2  °2  cos2  kx//dx) 

*  k  a/?/2  . 


If  we  set  (vn/n)r  „  s  -  a/L„.  l„  ■  If1.  ) 

Subject  to  this  uncertainty  the  various  graphical  results  of  Section  2  can 
be  Interpreted  In  terms  of  structure  scale  length  directly. 

The  results  from  Section  2  furnish  strong  limitations  on  the  persis¬ 
tence  to  lower  altitude  of  structure  due  to  fast  particle  structure  present 
at  the  bomb  burst  altitude.  Figures  11  and  13  establish  for  the  400  km  burst 
with  vertical  field  lines  that  degradation  Increases  with  decreasing  length 


t 
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scale,  also  that  the  major  contribution  to  degradation  comes  from  the  region 
In  altitude  where  the  neutral  mean  free  path  for  Impact  Ionization  Is  com¬ 
parable  to  the  structure  size.  At  higher  altitudes  charge  exchange,  propor¬ 
tional  to  the  ambient  neutral  density.  Is  relatively  weak;  at  lower  altitudes 
the  mean  free  path  for  Impact  Ionization  becomes  less  than  the  structure 
scale  length  and  diffusion  of  Ionized  structure  through  the  intervening 
neutral  transport  becomes  weaker.  The  degradation  of  structure  through 
neutral  transport  persists  to  lower  altitudes  for  shorter  wavelengths. 

Hence  the  Increased  damping  of  Fourier  components  with  decreasing  wave- 
1 ength. 

Figures  15  -  17  show  that  the  wavelength  dependence  of  degradation 
is  increased  when  the  magnetic  field  makes  an  angle  with  the  vertical.  This 
Is  due  to  the  greater  path  length  between  altitudes. 

The  results  of  Section  2.2  show  behavior  which  parallels  that  of 
Section  2.1.  However,  one  notes  that  Figure  24  Indicates  more  degradation 
than  corresponding  Figure  10.  This  appears  to  be  largely  due  to  the  greater 
densities  In  the  ionospheric  model  of  Section  2.2  which  result  In  the  region 
of  maximum  Increase  in  degradation  being  higher  in  altitude  where  the  atmos¬ 
pheric  density  gradient  scale  lengths  (and  hence  the  extent  of  region  of 
strong  degradation)  are  greater. 

The  figures  of  Section  2.3  show  for  fireball  behavior,  where  the 
cold  neutral  density  Is  severely  eroded,  that  degradation  from  neutral 
effects  is  less  severe. 

The  results  of  Section  3  show  that  chemical  recombination  affects 
structure  significantly  If  Hq  at  >  1.  (Here  ng  *  the  Ionization  at  the 
start  of  the  time  Interval,  t,  and  a  ■  the  recombination  coefficient.) 
Recombination  causes  all  structure  to  decay  more  rapidly  than  a  uniform 
background,  and  an  initial  structure  with  a  single  large  amplitude  Fourier 
component  produces  a  spectrum  of  higher  harmonics  which  exponentially  decay 
with  Increasing  wavenumber. 

The  results  of  Section  4  for  diffusion  due  to  Coulomb  scattering 
show  that  significant  degradation  from  this  mechanism  is  limited  to  wave¬ 
lengths  of  several  kilometers  or  less  provided  the  electron  density  Is  less 

Q 

than  10  cm  and  the  Ionization  temperature  Is  greater  than  0.1  eV. 
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The  results  of  this  report  support  the  potential  usefulness  of  the 
approach  of  separating  early  time  behavior  Into  two  Intervals: 

a)  0  <  t  <  tj,  when  plasma  wave  mechanisms  develop  structure. 

b)  tj  <  t  <  tg.  when  structure  decays  according  to  non-plasma 
wave  mechanisms,  prior  to  the  time  when  the  magnetic  field 
relaxes  back  to  ambient. 

It  Is  possible  due  to  the  strength  of  the  mechanisms  of  Sections  2  and  3 
that  detailed  consideration  of  much  of  the  plasma  Instability  behavior  which 
has  occurred  In  connection  with  the  early-time  behavior  of  HANE  (prior  to  t^) 
is  unnecessary  for  the  provision  of  Initial  condition  Inputs  to  the 
SCENARIO  code2  (at  t2). 

The  following  areas  and  questions  derive  from  this  report: 

1.  Detailed  application  of  these  notions  to  large  scale  codes. 

2.  Estimation  of  tj  as  a  function  of  free-energy  source  (When  do 
various  free  energy  sources  of  instability  turn  off?  Do  the 
turn  off  times  occur  prior  to  the  relaxation  of  the  magnetic 
field  back  to  ambient?). 

3.  Estimation  of  field  line  coupling  effects  (space-charge  and 
electromagnetic)  on  the  spectrum  of  structure. 

4.  Estimation  of  mass  flow  effects  on  structure  scale  lengths. 
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APPENDIX  A 

CLASSICAL  PLASMA  AND  MAGNETIC  FIELD  TRANSPORT 
IN  SLAB  GEOMETRY 


Two  mechanisms,  classical  and  anomalous,  relax  plasma  and  magnetic 

field  configurations  In  such  a  way  that  entropy  is  maximized  and  thermal 

A1-A4 

equilibrium  is  approached.  The  difference  between  the  two  mechanisms 

is  that  those  which  are  classical  are  a  result  of  particle  collisions  whereas 
those  which  are  anomalous  are  caused  by  wave  phenomena.  In  general,  the  mini¬ 
mum  rate  of  plasma  and  magnetic  field  transport  possible  in  any  configuration 
is  classical  with  appropriate  account  being  taken  of  geometry. 

In  this  appendix  the  rate  of  plasma  and  magnetic  field  transport  in 
a  slab  geometry  is  derived  assuming  classical  electron-ion  collisions.  The 
precise  plasma  model  used  In  the  calculations  is  shown  in  Figure  A-l.  The  z- 
dlrected  magnetic  field  of  strength  BQ  and  the  particle  pressure,  P,  have  x- 
directed  gradients.  The  y-dlrection  Is  ignorable. 

The  equations  which  give  the  plasma  and  magnetic  field  transport  In 
ti,'  x-direction  are  for  simply  Ionized  ions:^1"^ 


de  1  ■+ 

0  ■  TT  -  ?  7  • 


°"i  1  -*■ 

0  *  -J-  +  i.  7  j 

3t  e  v  Ji 


-*■  J’o  x  Bn  rn  +  ^ 

0  .  -vPe  -  <V  ♦  — T-*  ♦  f  (J,  ♦  J,) 

0  -  -VP,  +  e„,E  .  — -  —  v„  (J,  ♦  J,) 


0  -  7  *  -  T  <Je  *  !l> 


0  ■ 7  * E  sr  ■ 


(A-la) 


(A-lb) 


(A-lc) 


(A-ld) 


(A-le) 


(A-lf) 


In  Eqs.  (A-la-A-lf)  the  symbols  ng,  n^ ,  Pg,  P^,  E,  jg,  j^,  and  vg1  represent 
the  electron  density.  Ion  density,  electron  pressure,  electric  field,  ion 
pressure,  electron  current.  Ion  current,  and  electron-ion  collision  frequency, 
respectively.  The  constants  c,  mg  and  e  are  the  speed  of  light,  electron  mass 
and  the  magnitude  of  the  electron  charge,  respectively.  Also,  time  is  repre¬ 
sented  by  "t." 

Several  features  of  Eqs.  (A-la)  to  (A-lf)  are  noteworthy.  First,  the 
addition  of  Eqs.  (A-lc)  and  (A-ld)  and  using  Eq.  (A-le)  results  in  the  equili¬ 
brium  condition 

e(ni  -  ne)E  =  v(Pg  +  P.  +  BQ2/8ir)  .  (A-2) 

This  condition  reduces  to 

0  *  v(Pfi  +  P.  +  B02/8tt)  (A-3) 


if 


E  •  EQy  (A-4) 

is  a  y-directed  inductive  electric  field.  Second,  inertial  terms  have  been 
neglected  in  Eqs.  (A-lc)  and  (A-ld).  When  such  terms  are  significant,  non¬ 
linear  phenomena  (e.g.,  shocks)  rather  than  diffusive  transport  are  described 
A2 

by  the  equations. 

The  electron  and  Ion  currents  can  be  determined  from  Eqs.  (A-lc)  and 
(A-ld)  In  terms  of  the  small  parameter  ve-j / 1 ne I  where 

ne  =  "e  W  (A-5) 

is  the  electron- cyclotron  frequer  To  lowest  order. 


V  */• ;) 

(A-6a) 

O'  O' 

(A-6b) 
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where  x  Is  the  x-dlrected  unit  vector.  To  next  order. 


c  _a_ 

57  ax 
0 


(p,- 


+  Pe) 


(A-7) 


Note  that  the  convective  EQ  x  BQ  term  in  Eqs.  (A-6a)  and  (A-6b)  does  not 
explicitly  contribute  to  the  current  in  Eq.  (A-7).  Also,  Eq.  (A-7)  does  not 
depend  on  the  existence  of  an  ambipolar  electric  field. 

Substitution  of  Eq.  (A-7)  into  Eqs.  (A-la)  and  (A-lb)  leads  to  the 
following  expressions  for  the  rate  of  which  electron  and  ion  density  profiles 
change  with  time. 


an  an. 
_ e  _  _ i 

at  at 


» («0 


m_ 


'ei 


—  (p 

ax  v  e 


pi> 


_a_  n  _a_ 

ax  un  ax  e 


(A-8) 


with  the  diffusion  coefficient, 

Dn  ' (s;)  me  (Te  +  V  vei  •  <A‘9> 

In  deriving  the  second  form  of  Eq.  (A-8) 


p  ■  "e'Te  *  Ti> 

with  constant  electron  and  ion  temperature,  Tfi  and  ,  respectively. 

To  find  the  rate  at  which  the  magnetic  field  diffuses  with  time  it 
is  noteworthy  that  Eq.  (A-3)  implies  that  the  total  particle  plus  magnetic 
pressure  is  a  function  of  time  but  not  of  space;  i.e. , 

g(t)  -  Pe  +  Pi  +  B02/8tt  .  (A-10) 


The  differentiation  of  Eq.  (A-10)  with  respect  to  time  results  in 
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( A- 11 ) 
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o  _  dg(t) 
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From  Eqs.  (A-3),  (A-8),  and  (A- 11)  it  follows  that 

2  2 

JL  !°_  =  J_D  _L  !o_+  dg(t) 
at  8n  ax  n  ax  8n  dt 

or  since  the  function,  g(t)  does  not  depend  on  space 

where 


(A-12) 


(A-13) 


f ( x ,t)  =  Bq2/8u  -  g(t)  .  (A-14) 

Eqs.  (A-8),  (A-12)  and  (A-13)  indicate  that  the  diffusion  equations  for  the 
magnetic  pressure  and  density  have  identical  forms  only  if  g(t)  does  not  depend 
on  time  (i.e.,  the  total  pressure  in  the  system  is  a  constant  for  all  times). 
Furthermore,  the  directions  of  density  and  magnetic  field  transport  are  oppo¬ 
sitely  directed. 

To  determine  whether  density  or  magnetic  field  diffuse  more  rapidly 
it  is  useful  to  define  the  scale  lengths, 

Ln  s  |ne/(ane/ax) |  (A-15a) 

Lb  =  |BQ2/(aB02/ax)|  .  (A-15b) 

Equation  (A-3)  permits  the  scale  lengths  to  be  related: 


L 


n 


6LB 


( A- 16) 


with 


8  »  8tt  (Pe  +  P1  )/Bq2 


(A- 17) 


79 


The  scale  times  for  diffusion  of  density  and  f(x,t). 


Tn  =  lne/(3V3t^  (A-18) 

xf  =  |f(x,t)/[3f(x,t)/3t]|  ,  (A-19) 

are  then  related  from  iqs.  (A-8)  and  (A-13)  by 

xn/x^  a  0^  .  (A-20) 


For  g(t)  a  constant  Eq.  (A-20)  gives  the  ratio  of  the  time  scale  for  density 
variation  to  the  time  scale  for  magnetic  field  pressure  variation. 

It  follows  from  Eq.  (A-20)  that  the  density  diffuses  more  rapidly  than 
f(x,t)  for  e  <  1,  but  f(x,t)  diffuses  more  rapidly  than  the  density  for  0  >  1. 
Because  of  the  increasing  magnetic  field  modification  for  high  beta  plasmas, 
this  behavior  is  consistent  with  the  diamagnetic  properties  of  the  assumed 
plasma  model . 

In  addition  to  determining  the  time  and  spatial  evolution  of  the  density 
and  magnetic  field  it  is  possible  to  write  an  expression  for  the  evolution  of 
the  inductive  y-directed  electric  field.  From  Eqs.  (A-lf)  and  (A-12) 


±L  /_L  n  a  Bo  ,  dg(t)\ 
cBq  l  3x  n  3x  8tt  dt  J 


(A-21) 


and  so  the  electric  field  depends  not  only  on  the  magnetic  field  strength  and 
spatial  gradients  of  the  magnetic  field  but  also  on  the  time  variation  of  the 
total  pressure  within  the  system. 
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APPENDIX  B.  DECAY  OF  FAST  ION  STRUCTURE  IN  THE  IONOSPHERE 

ABSTRACT 

Fast  ions  with  energies  in  excess  of  ten  kilovolts  can  be  generated  with¬ 
in  the  ionosphere  by  magnetic  storms  or  manmade  mechanisms.  It  is  shown  that 
charge  exchange  of  fast  ions  into  fast  neutrals  results  in  the  rapid  attenu¬ 
ation  of  structure  with  neutral  mean  free  paths  much  larger  than  the  structure 
dimensions  and  the  generation  of  structure  of  larger  size.  For  structure  with 
dimensions  much  larger  than  the  neutral  mean  free  path  fast  neutral  loss  is 
diffusive  rather  than  convective  and  attenuation  of  structure  is  less  rapid. 


INTRODUCTION 

Fast  ions  with  energies  in  excess  of  one  kilovolt  can  be  generated  in  the 
ionosphere  by  magnetic  storms  [Sharp  et  al.t  1976a;  Sharp  et  al.,  1976b; 
Banks,  1979]  or  manmade  mechanisms  [Zinn  et  al.,  1966;  Glasstone  and  Dolan, 
1977].  Once  formed,  these  ions  travel  on  magnetic  drift  surfaces  toward  the 
northern  or  southern  conjugate  points.  While  following  their  trajectories  the 
ions  can  undergo  charge  exchange  collisions  with  less  energetic  background 
neutrals  [Fite,  et  al.,  1972]  thereby  generating  fast  neutrals,  with  energies 
approximating  the  original  fast  ion  energies,  and  slow  much  less  energetic 
ions.  The  fast  neutrals  are  not  constrained  to  magnetic  drift  surfaces  and 
follow  a  trajectory,  consistent  with  those  forces  which  are  not  electromagnet¬ 
ic  or  electrostatic  (e.g.,  gravity,  pressure,  etc.),  until  they  undergo  an  im¬ 
pact  ionization  collision  with  the  cold  neutral  background  [Fite,  et  al., 
1972].  With  the  impact  ionization  collision  a  fast  ion  and  a  less  energetic 
neutral  are  generated.  In  addition  to  undergoing  charge  exchange  and  impact 
ionization  processes,  the  fast  particles  can  concurrently  or  in  separate  col¬ 
lisions  ionize  cold  neutral  material  [Fite,  et  al.,  1972]  and  increase  the 
total  amount  of  ionization  available.  Successive  charge  exchange,  impact  ioni¬ 
zation  and  ionization  events  occur  until  the  energy  of  the  initially  fast  par¬ 
ticles  is  reduced  to  the  background  level. 
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The  purpose  of  this  appendix  is  to  calculate  the  attenuation,  resulting 
from  charge  exchange  and  impact  ionization  collisions,  of  an  initial  structure 
of  fast  ions  as  the  structure  travels  down  field  lines  into  denser  portions  of 
the  atmosphere.  It  is  demonstrated  that  the  attenuation  is  most  rapid  for 
small  structures  because  for  such  structures  the  neutral  mean  free  path  for 
impact  ionization  is  larger  than  the  structure  dimensions  and  neutral  trans¬ 
port  is  convective  rather  than  diffusive. 

Subsequent  portions  of  this  appendix  are  divided  in  the  following  way.  In 
the  second  and  third  sections  the  plasma  model  and  equation  are  described,  re¬ 
spectively.  In  the  fourth  and  fifth  sections  solutions  to  the  equations  in  a 
homogeneous  plasma  and  in  the  ionosphere  are  given.  The  final  section  is  a 
discussion  of  results. 


MODEL 

The  model  used  as  the  basis  for  the  calculations  is  shown  in  Figure  B-l. 
A  z-directed  magnetic  field  of  strength  BQ  is  assumed  antiparallel  to  the 
gravity  vector,  <J.  Hence,  the  magnetic  field  geometry  is  most  app  cable  to 
regions  near  the  north  or  south  magnetic  poles.  The  inhomogeneity  of  the  fast 
ion  density,  nj,  and  fast  neutral  density,  n2»  are  x-directed  and  the 
y-coordinate  is  ignorable.  The  evolution  of  an  infinitesimal  z-directed  length 
of  the  structure  is  followed  and  the  velocity  space  fast  ion  distribution 
function  is  assumed  for  definiteness  in  the  calculations  to  be  a  delta  func¬ 
tion  with  equal  perpendicular  and  parallel  speeds.  The  velocity  distribution 
of  neutrals  formed  by  charge  exchange  is  isotropic  in  the  x-  and  y-directions 
because  of  ion  gyration  about  the  magnetic  field.  To  determine  the  evolution 
of  a  finite  z-directed  length  of  structure  and  a  velocity  spread  of  particles 
solutions  analogous  to  the  ones  which  are  considered  here  would  have  to  be 
superimposed. 

The  fast  particles  are  taken  to  be  either  singly  ionized  or  neutral  mon¬ 
atomic  oxygen.  The  cold  neutral  background  species  and  density  as  a  function 
of  altitude,  n$,  are  taken  from  Fite  et  al.,  1972  .  Possible  variation  of  neu¬ 
tral  density  in  the  direction  orthogonal  to  the  magnetic  field  is  neglected. 


EQUATIONS  FOR  EVOLUTION  OF  FAST  PARTICLE  DENSITIES 

To  establish  the  evolution  of  the  fast  particle  structure,  the  calcula¬ 
tion  is  carried  out  in  a  Lagrangian  frame  of  reference  moving  with  the  initial 
fast  ion  structure  along  the  magnetic  field.  The  relevant  equations  are: 

dn, (x,t) 

- - -  -  v^tjn^x.t)  +  v2(t)n2(x,t)  (B-la) 

dn9(x,t)  H 

— ^ - =  -  v2(t)n2(x,t)  +  vl(t)n1(x,t)  -  v0(gy,t)  n2(x,t)  (B-lb) 

with  the  symbols  the  fast  ion  density  as  a  function  of  x  and  t  [nj(x,t)],  fast 
neutral  density  as  a  function  of  x  and  t  [n2(x,t)]  charge  exchange  collision 
frequency  as  a  function  of  time  [vj(t)3,  impact  ionization  collision  frequency 
[^(t)],  and  neutral  loss  frequency  which  is  a  function  of  a  differential  op¬ 
erator  and  time  [vQ(d/dx,t)].  It  is  useful  to  rewrite  Eqs.  (B-la)  and  (B-lb) 
in  terms  of  the  Fourier  integrals: 

oo 

nj(x,t)  =  /  dk  nj(k,t)  exp(ikx)  .  (B-2a) 

-00 

n2(x,t)  =  /  dk  n2(k,t)  exp(ikx)  .  (B-2b) 

Hence, 

+  Vj(t)  nj(k.t)  =  v2(t)n2(k,t)  (B-3a) 

$£+  v2(t)  +  v0(k,t)]n2(k,t)  =  v1(t)n1(k,t)  (B-3b) 
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The  collision  frequencies  In  Eqs.  (8-la),  (B-lb),  (B-3a)  and  (B-3b)  are 
determined  through  the  relations 


vl(t)  ■  V1  |  <’lsns't> 

(B-4a) 

v2(t)  =  v2  J  «2sns(t) 

(B-4b) 

vj,  V2»  ols,  o2S  and  n$(t)  represent  the  fast  ion  speed,  fast  neutral  speed, 
cross  section  for  collisions  of  fast  ions  with  cold  neutrals  of  species  s, 
cross  section  for  collisions  of  fast  neutrals  with  cold  neutrals  of  species  s, 
and  cold  neutral  density  of  species  s  as  a  function  of  time,  respectively. 

Two  forms  for  the  neutral  loss  frequency  have  been  chosen.  First  if  the 
neutral  mean  free  path  is  much  larger  than  the  structure  dimensions  in  the  x- 
direction,  neutrals  can  convect  out  of  the  structure  before  they  can  undergo 
an  impact  ionization  collision.  Hence,  such  particles  are  effectively  lost 
from  the  structure,  and  the  collision  loss  frequency  is  approximated  by  the 
following  expression  [Sears,  1953]. 

vQ(k,t)  =  kv2/8  (B-5) 

For  cases  when  the  neutral  mean  free  path  is  much  smaller  than  the  structure 
dimension  in  the  x-coordinate  a  diffusive  loss  model  is  appropriate  [Sears, 
1953]  and 


,  (I-  t)  -  k  X 
vo(k,t)  3ul 


(B-6) 


In  the  numerical  solution  of  Eqs.  (B-3a)  and  (B-3b)  the  expression  for  the 
collisional  loss  term  is  taken  to  be  the  minimum  of  Eqs.  (B-5)  and  (B-6). 
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SOLUTION  OF  EQUATIONS  FOR  A  HOMOGENEOUS  IONOSPHERE 


Qualitative  features  of  the  fast  particle  structure  decay  can  be  discov¬ 
ered  by  solving  Eqs.  (B-3a)  and  (B-3b)  in  the  limit  when  the  collision  fre¬ 
quencies  and  the  neutral  loss  frequency  are  time  independent.  Hence,  in  this 
section  Eqs.  (B-3a)  and  (B-3b)  are  solved  and  discussed  in  this  limit. 

Eqs.  (B-3a)  and  (B-3b)  can  now  be  written  as: 

[p(k)  +  vj ]n1(k)  =  v2n2(k)  (B-7a) 


[  P(  k)  +  vQ  +  v2]n2(k)  =  vjn^k) 


with 


nj(k,t)  ~  exp[p(k)t] 


(B-7b) 


(B-8a) 


n2(k,t)  ~  exp[p(k)t] 


(B-8b) 


It  follows  from  Eqs.  (8-7a)  and  (B-7b)  that  the  explicit  values  for  p(k)  are  a 
solution  of  the  characteristic  equation: 

0  =  p(k)2  +  (vj  +  v2  +  vQ)p(k)  +  vjv0  .  (B-9) 

Hence: 

P±(k) 

where  the  two  solutions  to  Eq.  (B-9)  have  been  written  as  p±(k)  depending  on 
the  sign  in  front  of  the  square  root  in  Eq.  (B-10).  Both  p±(k)  are  negative 
definite  and  contribute  to  decay  of  structure.  Based  on  Eqs.  (B-8a),  (B-8b) 
and  (B-10)  the  explicit  solution  to  Eqs.  (B-7a)  and  (R-7b)  can  be  written  as 


(V1  +  v2  +  v0)  *  ^V1  +  v2  +  ^  "  4vlvo^1/2 


(B-10) 
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n j (k,t)  = 


1 


pfk)  l  v'  A(k)  exP(P+t)  +  p~~T  Vl  B00  e*P(P.t)  ,  (3-lla) 

n2(k,t)  =  A(k)  exp(p+t)  +  B(k)  exp(p_t)  .  (B-llb) 

where  A(k)  and  B(k)  are  constants.  With  the  initial  conditions 

n,(k,t) I  =  n  (k)  (B-12a) 

1  't=0 

n?(k,t)|  =  0,  ( B- 12b ) 

*t=0 

Equations  (B-lla)  and  (B-llb)  can  be  rewritten  as 

n0(k) 

ni(kjt)  =  __ — -  {(p_  +  exp(p+t)  -  (p+  +  Vj)  exp(p_t)}  (B-13a) 
,  4  no00  (P+  +  v,)(p_  +  v.) 

n2(k,t)  »  p  : . - — -  [exp(p+t]  -  exp(p_t)]  (3-13b) 

with  the  dependence  of  p+  and  p_  on  wavenumber  suppressed  for  brevity  of  nota¬ 
tion. 

It  is  useful  to  rewrite  Eqs.  (8-13a)  and  (B-13b)  in  two  limits.  First, 
for  vj,  \>2  <  vQ  fast  neutrals  are  lost  almost  immediately  after  being  formed 
and  the  equations  become 

n^k.t)  -  nQ(k)  expf-ujt)  (3-14a) 

n2(k,t)  *  (vj/Vq)  nQ(k)  expf-Vjt)  .  (B-14b) 

An  important  feature  of  these  equations  is  that  fast  ion  structure  only  de¬ 
pends  on  its  initial  value  and  the  charge  exchange  collision  frequency.  Sec¬ 
ond,  for  >  vQ  the  total  number  of  fast  ions  plus  fast  neutrals  decreases 

very  slowly  relative  to  the  charge  exchange  and  impact  ionization  collision 
times  and 


n^k.t) 


(B- 15a ) 


n0(k) 

"  v'  V  v2  ^2  +  V1  exPHvl  +  v?)ti} 

n0(k) 

n2(k.t)  “  v"°V  v  V1  ^  •  exP[-Cvi  +  v2 )t ]  1  (3-  15b) 

Hence,  for  time  approaching  infinity 

n1(k,t)/n2(k,t)  -  v?/vj  .  (B-16) 

which  is  independent  of  k. 

SOLUTION  OF  EQUATIONS  FOR  AN  INHOMOGENEOUS  IONOSPHERE 

Equations  (B-3a)  and  (B-3b)  have  been  solved  numerically  for  the  real  in¬ 
homogeneous  ionosphere.  The  initial  conditions  for  the  calculation  are  an 
altitude  of  800  kilometers  at  t  =  0  and 


n0(k)  =  n. (k,t) 

1  t=0 

(B- 17a) 

0  =  n,(k,t) 

c  t=0 

(B-17b) 

The  numerical  solution  is  followed  down  the  field  lines  to  an  altitude  at 
which  the  fast  ion  density  for  each  wavelength  is  about  0.001  of  its  initial 
value.  For  vertical  magnetic  field  lines,  a  fast  particle  energy  of  12  keV, 
and  equal  parallel  and  perpendicular  speeds,  the  fast  ion  velocities  down  the 
field  line  are  1.9  x  107  cm/sec.  The  calculations  have  been  carried,  out  for 
two  different  ionospheric  models  with  the  cold  neutral  diatomic  nitrogen,  di¬ 
atomic  oxygen  and  monatomic  oxygen  densities  at  each  altitude  given  by  the 
C IRA  1965  Hour  10,  Model  9  and  1965  Hour  22,  Model  4  tables,  respectively  [CO¬ 
SPAR  Working  Group  IV,  1965].  (The  first  model  is  appropriate  to  a  quiescent 
ionosphere  with  a  very  high  level  of  solar  activity  while  the  second  model  is 
appropriate  to  an  ionosphere  with  a  level  of  mean  solar  activity  but  under  the 
magnetic  storm  conditions  of  Sharp,  et  a!.,  1976a  and  Sharp,  et  al.,  1976b.) 
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Results  are  summarized  in  Figures  B-2  and  B-3  for  wavenumbers  with  k  > 
10”^  cm"*  and  the  two  ionospheric  models.  The  calculations  indicate  that  the 
fast  ion  attenuation  is  insensitive  to  the  wavenumber  particularly  at  higher 
altitudes.  Ion  structure  decays  by  three  orders  of  magnitude  in  going  from  8ft0 
km  to  370  km  in  altitude  in  Figure  B-2  and  from  800  to  280  km  in  Figure  B-3. 
Also,  the  fast  neutral  density  is  calculated  to  be  appreciably  smaller  than 
the  fast  ion  density.  These  results  are  consistent  with  the  homogeneous  cold 
neutral  results  of  the  previous  section  in  the  limit  when  the  mean  free  path 
for  neutral  collision  is  much  larger  than  structure  dimensions  [Eqs.  (B-14a) 
and  (B-14b)].  Hence,  to  lowest  order  the  structure  behavior  is  more  closely 
approximated  by  Eqs.  (B-14a)  and  (B-I4b)  rather  than  Eq.  (B-16)  and  the  evolu¬ 
tion  of  the  fast  ion  structure  in  the  Lagrangian  frame  of  reference  moving 
along  the  magnetic  field  with  the  fast  ions  is 

t 

n. (k,t)  =  n  (k,t)  exp[-  /  dt  v,(t)]  .  (B-18) 
io  0 

For  lower  altitudes  or  larger  wavelengths  the  neutral  mean  free  path  is 
much  smaller  than  the  structure  size  and  the  loss  of  particles  is  governed  by 
neutral  diffusion  rather  than  neutral  convection.  It  follows  from  Eq.  (B-6) 
that  this  diffusive  loss  is  weakest  for  inverse  wavenumber  and  impact  ioniza¬ 
tion  frequency  going  to  infinity.  Under  circumstances  when  v0(t)  <  v^(t), 
^(t)  the  relationship  between  nj(k,t)  and  n2(k,t)  is  governed  by  an  adiabatic 
relation  analogous  to  Eq.  (B-18), 


n1(k,t)/n2(k,t)  =  v^t)/^)  .  (B-19) 


In  the  limit  of  k  0,  the  relation  between  the  fast  ion  and  fast  neutral  den¬ 
sities  as  a  function  of  the  initial  ion  density,  Eq.  (B- 17a) »  is: 


n0(k) 

ni(k,t)  =  1  +  vjWv'gTTT 


(B-20a) 


n2(k,t) 


[v1(t)/v2(t)]nQ(k) 
1  +  v1(t)/v2(t) 


(B-20b) 
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Eq.  (B-20a)  indicates  that  at  the  altitude  where  the  ion  density  for 
k  +  •  is  one  thousand  times  smaller  than  its  initial  value  (i.e.,  370  km  for 
Figure  B-2  and  280  km  for  Figure  B-3)  the  ratio 


n j (k,t) 

k+» 

nj(k,t) 

k+0 

10  ^  [1  + 


(B-2 1 ) 


For  the  neutral  densities  and  cross-sections  appropriate  to  both  Figures  B-2 
and  B-3  this  ratio  is  approximately  7  x  10"^. 

It  is  clear  that  the  exponential  decay  rate  for  structure  at  a  given 
wavenumber  is  greatest  for  altitudes  where  the  neutral  mean  free  path  is  com¬ 
parable  to  the  inverse  wavenumber. 


DISCUSSION 

Fast  ions  can  be  generated  by  both  geophysical  or  manmade  phenomena  in 
the  upper  atmosphere. 

It  has  been  shown  that  the  charge  exchange  of  fast  ions  (e.g.,  monatomic 
oxygen  with  12  keV  kinetic  energy)  into  fast  neutrals  is  a  powerful  mechanism 
for  attenuating  fast  ion  structure  from  the  altitude  of  8fio  km  down  to  370 
kilometers  (for  Figure  B-2)  and  280  kilometers  (for  Figure  B-3)  as  long  as 
Eqs.  (B-14a)  and  (B-14b)  are  a  good  description  for  the  decay.  For  the  parame¬ 
ters  and  ionospheric  model  appropriate  to  Figures  B-2  and  B-3  this  corresponds 
to  an  inverse  wavenumber  of  approximately  2  and  1  km,  respectively.  For  larger 
inverse  wavenumbers  or  lower  altitudes  diffusion  of  fast  neutrals  becomes  more 
important  and  so  the  wavenumber  component  of  structure  decays  at  a  less  rapid 
rate.  Consequently  for  structure  consisting  of  significant  wavenumber  compo¬ 
nents  with  inverse  wavenumbers  less  than  or  greater  than  the  ion  and  neutral 
mean  free  path,  decay  as  well  as  structure  shape  changes  occur  due  to  strong 
attenuation  of  the  short  wavelength  contribution. 

The  calculations  also  demonstrate  that  the  penetration  of  fast  ion  struc¬ 
ture  originating  at  high  altitude  to  lower  altitudes  is  a  strong  function  of 
ionospheric  model.  The  physical  basis  for  this  effect  is  the  increase  in  neu¬ 
tral  densities  especially  at  altitudes  above  300  km  with  increasing  solar 
activity. 


9 
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Figure  B-2.  Fast  ion  density  for  k  -*■  »  in  units  of  the  initial  fast 
ion  density  at  800  km  vs.  altitude  in  kilometers. 
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APPENDIX  C.  ELECTRIC  FIELD  LIMITATION  OF  DRIFT  DISSIPATIVE 
WAVE  GENERATION  IN  IONOSPHERIC  STRUCTURES 

ABSTRACT 

Coupled  background  electric  fields  and  density  profiles  exist  across  the 
sides  of  striations  and  other  elongated  ionospheric  structures.  The  presence 
of  the  electric  field  profoundly  modifies  the  behavior  of  drift  dissipative 
waves  by  shifting  the  location  of  the  waves  away  from  the  point  where  the  rel¬ 
ative  density  gradient  is  a  maximum,  consequently  inhibiting  the  growth  of  the 
waves. 


INTRODUCTION 

Plasma  waves  in  the  ionosphere  are  sources  of  electron  density  structure. 
In  addition  within  the  laboratory  it  is  well  known  that  they  can  lead  to  anom¬ 
alous  diffusion  of  pre-existing  density  structures  [Cheng  and  Okuda,  1977]. 

It  has  been  observed  that  barium  cloud  strations  and  other  ionospheric 
plasma  density  structures  (as  in  Equatorial  Spread  F)  tend  to  be  elongated  in 
the  plane  perpendicular  to  the  geomagnetic  field  in  the  direction  of  the  ?  x  § 
drift  [Perkins  and  Doles,  1976;  Scannapieco  and  Ossakow,  1976;  Tsunoda  et  al., 
1979].  (Here  P  and  P  are  the  ambient  electric  and  magnetic  fields,  respective¬ 
ly.)  Around  most  of  the  perimeter  of  such  structures,  the  ambient  electric 
field  should  be  parallel  (or  anti-parallel)  to  the  density  gradient.  Gradient- 
drift  modes  have  been  shown  to  be  stabilized  under  these  conditions  [Perkins 
and  Doles,  1976],  Hudson  and  Kennel  [1975]  have  pointed  out  the  possible  Im¬ 
portance  of  drift-dissipative  modes  in  ionospheric  structures.  The  purpose  of 
this  calculation  is  to  evaluate  the  effect  of  this  electric  field  on  drift 
dissipative  mode  growth  rates.  The  key  result  of  the  calculation  is  that  drift 
dissipative  modes  tend  to  be  stabilized  through  the  presence  of  ambient  elec¬ 
tric  fields. 
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This  appendix  Is  divided  Into  two  additional  sections.  In  the  second  sec* 
tion,  the  dispersion  relation  for  drift  dissipative  modes  In  the  presence  of 
an  ambient  electric  field  Is  derived  and  solved.  The  final  section  Is  a  short 
discussion  and  review  of  the  results. 


DISPERSION  RELATION 


a.  General  Analysis 

In  this  section,  the  nonlocal  dispersion  relation  for  drift  dissipative 
modes  In  the  presence  of  an  ambient  electric  field  is  derived  and  solved.  As 
the  basis  for  the  calculation,  a  plasma  model  Is  used  (see  Figure  C-l)  In 
which  the  x-directlon  is  normal  to  the  surface  of  the  striation  perimeter.  The 
plasma  density  gradient  and  the  ambient  electric  field  are  then  x-directed. 
The  plasma  Is  assumed  to  be  homogeneous  in  the  y-  and  z-di recti ons  and  a  uni¬ 
form  z-dlrected  field  is  present.  The  plasma  itself  consists  of  a  quasi¬ 
neutral  mixture  of  electrons  and  singly  ionized  ions.  A  uniform  background  of 
neutral  particles  is  also  permitted  in  the  model. 

Consistent  with  the  plasma  model,  the  wave  electric  field  has  the  form 


^(x.y.z.t)  =  -v{*j(x)  exp[1(kyy  +  k2z  -  ut)]}  ,  (C-l) 


where  t  denotes  the  time.  In  a  colllsionless  plasma  and  in  the  absence  of  the 
ambient  electric  field,  the  dispersion  relation  for  drift  waves  with  frequen¬ 
cies  well  below  the  ion-gyrofrequency  has  the  form  [Cheng  and  Okuda,  1976] 


A, 


0  , 


(C-2) 


where 


> 


o  -  *'*  *,[1  +  <•  *  ■tf’W71.  • 

4  *  Pf (»  +  -  Xj^)/Tj  t 

<  ■  -ky  *k  «•[!»(»)  ]«*  ■ 

*«d  ■  [,<ba>  «P(-ba)  *(<■/*  kzva)/^k2va  , 
ba  *  V/  * 

Va  ’  (T«/m<,)1/2  • 

fl  3  q  B„/m  c  . 

a  a  0  o 


(C-3) 


The  magnetic  field  strength,  the  speed  of  light,  the  charge  of  a  species  a 
particle,  the  temperature  of  species  a,  and  the  mass  of  species  a  are  here 
represented  by  BQ,  c,  qa,  Ta,  and  ma,  respectively.  Electron  and  ion  parame¬ 
ters  are  differentiated  by  subscripts  e  and  i.  The  plasma  dispersion  function 
(Fried  and  Conte,  1961)  and  the  modified  Bessel  function  of  the  first  kind  are 
denoted  by  Z  and  It,  respectively.  Also,  the  dependence  of  Q,  A,  and  n  on 
x  has  been  suppressed  in  the  notation  of  Eqs.  (C-2)  and  (C-3). 

To  modify  Eq.  (C-2)  into  a  form  appropriate  to  drift  dissipative  modes  in 
the  presence  of  an  ambient  electric  field,  the  following  replacements  are  made 
[Kadomtsev,  1965]: 

u'  =  u  +  kycE0(x)/B0  ,  (C-4a) 


[1  +  (“  +  “e)X0eJ/Te  *  TJ  (  +  f«)  • 

where 

«  *  (kzve)2/ve  ‘ 


(C-4b) 

(C-5) 
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In  Eqs.  (C-4)  and  (C-5),  the  quantities  E0(x)  and  ve  denote  the  ambient  elec¬ 
tric  field  strength  and  the  electron  collision  frequency,  respectively.  For 
electron  collisions  with  ions, 

23.4  -  1.15  log1Q(n  )  +  3.45  l°9io(Te)  n 

ve  * - - 5 - -T77  .  (C-6) 

3.5  x  10b  T3/Z 

e 

where  Te  and  ne  are  in  units  of  electron  volts  and  inverse  cubic  centimeters, 
respectively. 

To  Illustrate  the  effect  of  the  electric  field  on  the  drift  mode,  it  is 
assumed  that 


«»M.  I«*|,  |w*|  .  (C-7a) 

In  addition, 

H  >  (2)1/2kzvi  (C— 7b) 

for  drift  waves,  and  so  the  parameter  Q  in  Eq.  (C-3)  can  be  written  as 

q(x)  -y-i-  [aw'  -  but]  ,  (C-8a) 

where 

a  *  1  +  Ti/Te  -  IqO^)  exp(-b^)  ,  (C-8b) 

b  -  IQ(b1 )  exp(-b^)  .  (C-8c) 


Now  If  the  x-directed  scale  length  of  the  wave  is  much  smaller  than  the 
density  gradient  scale  length. 


(C-9) 


then  the  parameter  Q  can  be  expanded  In  a  Taylor  series  about  the  point,  x«0, 
where  the  local  dispersion  relation  is  approximately  satisfied;  i.e.. 


d2q^x) 


(C-10) 


x-0 
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The  first  term  on  the  right-hand  side  of  Eq.  (C-10)  Is  the  local  dispersion 
relation  for  drift  waves  when  equated  to  zero  [Kadomtsev,  1965].  The  last  two 
terms  contribute  to  the  nonlocal  dispersion  of  the  modes. 

Now  If 

Q(x) |x=0  *  0  ,  (C-lla) 


x*0 


0  , 


then 


d2Q(«> 

dX 


1  -  dw'  l. 

„0  ■T7=-4ar-  b3x- 


x-0 


x»0 


i 

f1A“  d?  dir 


x*0 


(C-llb) 


(C- 12a) 

(C-12b) 


In  Sub-appendix  CA,  relations  [Eqs.  (C-8a)-(C-9b)]  between  the  background 
density  gradient  and  electric  field  strength  are  evaluated.  Using  these  ex¬ 
pressions,  Eqs.  (C-12a)  and  (C-12b)  can  be  rewritten  In  the  following  way. 


(C-13b) 


b.  Asymptotic  Analysis 

Let  denote  the  maximum  value  of  u?  [and  hence  of  (l/n)(dn/dx)]  over 

1  flH  » 

the  range  In  x  with  spatial  variation. 
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The  effect  of  the  background  electric  field  can  be  categorized  by 


r 


(C-14) 


For  r  >  1  electric  field  effects  In  Eqs.  (C-13a)  and  (C-13b)  are  Insignifi¬ 
cant,  the  drift  mode  is  localized  around  the  value  of  x  such  that  |«*  |  Is  a 
maximum,  and  the  growth  rate  scales  as  |.  (See  Sub-appendix  CB  *for  de¬ 
tails.) 


For  r  <  1  the  electric  field  effect  Is  dominant.  This  limit  obtains  for 
density  gradient  scale  lengths  perpendicular  to  B  (■  Ln)  with 


(C-15) 


Typically  v^/(cEX0/B)  Is  less  than  10  and  the  inequality  is  satisfied  for 
Ln  >  0.1  km. 

To  demonstrate  the  electric  field  effect  quantitatively  for  r  <  1,  assume 

that 


i  flfn 

n  dx< 


(C- 16) 


with  £  of  order  unity,  and  the  maximum  value  of  [l/n(dn/dx)]  denoted  by 
[l/n(dn/dx)]m. 

Condition  (C-16)  appears  to  be  likely  for  higher  order  density  deriva¬ 
tives,  although  one  could  mathematically  construct  density  profiles  which 
violate  it.  On  using  Eqs.  (C-llb)  and  (C-13a)  one  obtains: 


vRi  i1  *  ¥Te) 

aky(ce-/BoJ(n-/nJ 
7*  / 1  dn  \ 

,  Vlac1\n  •ax/x-0 
"  a  ky(cEtB/Bj(nt./n) 


[l  ♦  (Ti/Tfi)] 


(C— 17) 


with  a  of  order  unity  for  T^  »  Te. 
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Using  Eq.  (C-16)  to  lowest  order  In  r  this  becomes: 


From  Eqs.  (C-16)  and  (C-18)  It  Is  clear  that 


(C-18) 


(C- 19) 


c.  Nonlocal  Behavior 

We  note  additionally  that  previous  nonlocal  treatments  of  drift  waves 
have  been  In  the  zero  ambient  electric  field  limit  and  have  been  carried  out 
about  the  points  where  the  magnitude  of  the  diamagnetic  drift  frequency  is  a 
maximum  and  d2q/dx2  Is  positive.  However,  with  an  ambient  electric  field,  Eq. 
(C— 19)  Indicates  that  the  drift  waves  are  centered  about  points  which  are  dif¬ 
ferent  from  the  ones  at  which  the  magnitude  of  the  drift  frequency  is  a  maxi¬ 
mum.  Moreover,  the  extremum  values  of  Q  are  shifted  in  space  as  a  result  of 
the  ambient  electric  field  and  It  Is  not  obvious  that  d2q/dx2  is  positive  or 
negative.  Hence,  in  Sub-appendix  CC  the  nonlocal  dispersion  of  drift  waves  is 
examined  first  for  d2q/dx2  greater  than  zero  and  then  less  than  zero. 

It  should  be  noted  that  the  detailed  features  of  the  x-directed  mode 
structure  are  higher  order  in  the  inverse  density  gradient  scale  length  than 
the  shift  in  modal  location.  Consequently,  these  detailed  features  are  corre¬ 
spondingly  less  significant. 


DISCUSSION 

A  nonlocal  analysis  of  the  dispersion  relation  for  drift  dissipative 
modes  in  the  presence  of  an  ambient  electric  field  has  been  evaluated.  A  crit¬ 
ical  feature  of  the  calculation  Is  the  coupling  of  the  ambient  electric  field 
to  the  density  profile  (Sub-appendix  CA). 
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The  analysis  suggests  that  drift  dissipative  waves  are  In  general  stabil¬ 
ized  by  the  presence  of  an  ambient  electric  field  which  shifts  the  location  of 
the  drift  waves  away  from  the  points  where  the  relative  density  gradient 
(l.e.,  n"/n)  Is  a  maximum  and  the  growth  of  the  waves  Is  most  vigorous.  In 
particular,  for  situations  where  the  E  x  B  velocity  Is  larger  than  the  diamag¬ 
netic  velocity  at  the  location  of  maximum  logarithmic  density  gradient  by  a 
factor  of  r  (»1),  typically  for  gradient  scale  lengths  Ln  >0.1  km,  the 
growth  rate  for  instability  Is  a  factor  of  r“*  lower  than  the  expected  value 
In  the  absence  of  the  ambient  electric  field. 

Further  these  calculations  suggest  that  drift-dissipative  waves  could  be 
triggered  on  the  sides  of  elongated  F  region  structures  by  the  increase  of  at¬ 
mosphere  ionization  in  the  E  region  at  dawn.  The  creation  of  a  highly  conduct¬ 
ing,  relatively  uniform  (iB)  region  of  ionization  connected  by  magnetic  field 
lines  to  the  altitude  of  structure  could  act  to  short  out  background  electric 
fields  acting  along  the  sides  of  structure.  The  resulting  growth  rates  for 
drift-dissipative  waves  would  then  be  substantially  increased.  This  could  lead 
to  anomalous  diffusion  of  Spread  F  structure  in  a  manner  analogous  to  that  al¬ 
ready  discussed  in  connection  with  laboratory  plasmas  [Cheng  and  Okuda,  1977]. 


SUB-APPENDIX  CA 


In  this  sub-appendix,  a  relationship  between  the  ambient  electric  field 
and  density  profile  is  derived  based  upon  the  constraint  that  the  flow  across 
the  density  gradient  is  quasi-neutral. 


Neglecting  terms  which  do  not  contribute  to  the  formation  of  amblpolar 
electric  fields  and  Inertial  terms,  the  Ion  momentum  equation  has  the  form 


Vi  * 

m^c 


-  Vl  ♦ 


m. 


m^n 


vn 


(CA- 1) 


where  vj  Is  the  ion-neutral  collision  frequency. 


1M 


In  an  expansion  In  the  small  parameter  Vj/Oj,  the  lowest-order  component 
of  the  ion  velocity  In  the  direction  of  the  density  gradient  Is 


v1x<x> 


In  the  Ionosphere, 


qEoxW  Ti  dn(x)  V1 
m^  "  m^n(x)  dx  "^7  * 

|ve/oel  <  |  , 


(CA-2) 


(CA-3) 


where  ve  is  the  electron- neutral  collision  frequency.  Hence,  the  quasi-neutral 
condition  requires  that 

const.  *  n(x)  v^x(x)  .  (CA-4) 


For 


0  *  1imfdn(x)/dx]  , 

X-H» 


(CA-5) 


and  assuming  that  only  n(x)  and  Eox(x)  are  functions  of  position,  then 

«Eox<x>  1  dn(x)  _  qE»n» 

m^  “  mi  h(xT  dx  ~  n(x)mi  * 


(CA-6) 


where 

E  ■  llmfE  (x)]  ,  (CA-7a) 

•  x-**»  ox 

nw«J1»[n(x)]  .  (CA-7b) 

It  follows  from  Eq.  (CA-6)  thet 


cE  n 


♦  df(»)  •  •  ♦  ky  . 

(CA-8) 

and 

dw'(x) 

(CA-9a) 

dx 

.  A*  eE.".  df  1  Mill 

(CA-9b) 
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SUB-APPENDIX  CP 


It  has  been  shown  that  a  large  ambient  x-directed  electric  fie'd  shifts 
the  location  of  drift  wave  localization  away  from  the  location  of  largest  den¬ 
sity  gradient.  To  show  that  the  shift  causes  a  reduction  of  drift  wave  growth 
rate  from  its  optimum  value,  the  local  dispersion  relation  with  the  destabil¬ 
izing  parallel  electron  diffusion  term  can  be  written  in  the  following  way: 


where 


0  * 


li  1  B(a  -1)  +  f»  -  g2) 
Te  a  +  g f 

a  =  -u'/u*  ,  6  = 


+  a  - 


-6/w£ 


(CB-1) 


(CB-2) 


Fixing  the  value  of  ky  and  permitting  k|  to  vary  in  such  a  way  that  0  is  con¬ 
stant,  it  is  evident  that  the  solution  for  the  real  and  imaginary  value  for 
scales  like  (i.e. , 

w"  ~  a»*d  (CB-3) 


where  d  is  a  constant).  Hence,  there  is  a  growth  rate  reduction  by  the  ratio 
of 

*  —I —  (CB-4) 


e 

RT 


max 


'max 


at  a  location  where  the  drift  frequency  is  reduced  from  its  maximum  value. 
Here,  the  quantity  y/vmx  is  the  ratio  of  the  drift  wave  growth  rates  at  the 
location  of  drift  wave  localization  with  and  without  an  electric  field,  re¬ 
spectively. 


In  the  preceding  discussion  the  stabilizing  role  of  ion-neutral  colli¬ 
sions  has  not  been  considered.  However,  Kadomtsev  points  out  that  for  ion- 
neutral  collision  frequency,  v^. 


|  <  1 


(CB-5) 
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leads  to  wave  stabilization  [Kadomtsev,  1965].  Hence,  for  fixed  ion-neutral 
collision  frequency,  reduction  of  drift  wave  frequency  from  its  maximum  value 
with  fixed  ky  leads  to  a  reduction  of  parameter  space  permitting  drift  wave 
growth.  It  follows  that  a  large  ambient  electric  field  lessens  the  vigor  of 
the  drift  waves  both  with  and  without  stabilization  by  ion-neutral  collisions. 


SUB- APPENDIX  CC 


In  this  sub-appendix  Eq.  (CC-2)  is  explicitly  solved  for  d2Q/dx2  greater 
than  and  less  than  zero,  respectively. 


1.  d2q/dx2  >  0 

Through  the  use  of  Eqs.  (CC-10),  (CC-llb)  and  (CC- 13b) ,  Eq.  (CC-2)  can  be 
written,  for  d2Q/dx2  >0,  in  the  form 


0  = 


1  .  _  u 

7+  P  “  J~ 


♦i(u) 


where 


u  =  x  exp(-iir/4)/e 


(CC-1) 


(CC-2) 


Q 


XsQ  =  -i ( 1/2  +  p)  B-2 


(CC-3) 


=  -1/2  B'4  • 
x=0 


(CC-4) 


The  appropriate  boundary  conditions  for  the  differential  equation  correspond 
to  the  outward  propagation  of  energy  away  from  x  =  0.  Hence,  the  index  p  has  a 
non-negative  integral  value,  and  the  solution  to  Eq.  (CC-1)  is  the  parabolic 
cylinder  function  [Morse  and  Feshbach,  1953], 


♦l(x)  =<  Dp[x  exp(-iir/4)/B]  • 


(CC-5) 


From  Eqs.  (CC-3)  and  (CC-4)  It  follows  that  there  Is  a  negative  contribution 
to  the  growth  of  the  wave  caused  by  the  outward  convection  of  energy  away  from 
x  -  0  CPearl stein  and  Berk,  1969;  Sperling  and  Perkins,  1976].  The  explicit 
form  for  this  negative  contribution  is 


Yd  =  -I*»*l(l  +  2p)|l  +  -£■  j[IQ(b.j)  -  (bi )  ]  exP('bi) 

r3  )1/2 

xexp(-b-)  [2o|  -  3013]  j 
2.  d2Q/dx2  <  0 

For  the  case  that  d2Q/dx2  <  0,  Eq.  (CC-2)  becomes 

°  “[d7  +  *  +  p  "2T*l(u)  ’ 

where 

u  =  x/8  , 

Q|x.0  -  (1/2  +  P)B'2  , 


1  +  ”  ^g(bi  ) 


(CC-6) 


(CC-7) 

(CC-8) 

(CC-9) 


d2Q 

dir 


x=0 


■1/2  e 


-4 


(CC-10) 


The  appropriate  boundary  conditions  for  Eq.  (CC-7)  is  that  ^  is  localized 
around  the  point  x  *  0  [l.e. ,  $j(u)  +  0  for  |x|  •].  The  solution  to  Eq. 

(CC-7)  Is  the  parabolic  cylinder  function 

4j(x)  *  Dp(x/&)  (CC-11) 

for  p  equal  to  a  non-negative  Integer. 
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From  Eqs.  (CC-2)  and  (CC-3),  It  follows  that  there  Is  a  real  frequency 
shift,  6u,  associated  with  the  nonlocal  dispersion  relation: 

■  -u* | Yd/»*|  (CC-12) 

where  Is  given  by  Eq.  (CC-6). 
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